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Abstract 

Energy for a nucleus is considered in macroscopic limit, in terms of nucleon numbers. Further 
considered for a nuclear system is the Hohenberg-Kohn energy functional, in terms of proton and 
neutron densities. Finally, Skyrme-Hartree-Fock calculations are carried out for a half- infinite 
particle-stable nuclear-matter. In each case, the attention is focused on the role of neutron-proton 
asymmetry and on the nuclear symmetry energy. We extend the considerations on the symmetry 
term from an energy formula to the respective term in the Hohenberg-Kohn functional. We show, in 
particular, that in the limit of an analytic functional, and subject to possible Coulomb corrections, 
it is possible to construct isoscalar and isovector densities out of the proton and neutron densities, 
that retain a universal relation to each other, approximately independent of asymmetry. In the so- 
called local approximation, the isovector density is inversely proportional to the symmetry energy 
in uniform matter at the local isoscalar density. Generalized symmetry coefficient of a nuclear 
system is related, in the analytic limit of a functional, to an integral of the isovector density. 
We test the relations, inferred from the Hohenberg-Kohn functional, in the Skyrme-Hartree-Fock 
calculations of half-infinite matter. Within the calculations, we obtain surface symmetry coefficients 
and parameters characterizing the densities, for the majority of Skyrme parameterizations proposed 
in the literature. The volume-to-surface symmetry-coefficient ratio and the displacement of nuclear 
isovector relative to isoscalar surfaces both strongly increase as the slope of symmetry energy in 
the vicinity of normal density increases. 

PACS numbers: 21.10.Dr, 21.10.Gv, 21.60..Jz, 21.65.-f, 21.65. Cd, 21.65.Ef 

Keywords: symmetry energy, half-infinite matter, nuclear matter, Hohenberg-Kohn functional, Skyrme- 
Hartree-Fock model, nuclear surface, isovector density, surface symmetry coefficient 
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I. INTRODUCTION 



Nuclear symmetry energy ties different areas of nuclear physics, including structure of 
ground-state nuclei, dynamics of central nuclear reactions, physics of giant collective exci- 
tations and physics of neutron stars. Recently, the interest in symmetry energy has been 
stirred up by novel astrophysical observations and by the availability of exotic beams in 
accelerators, with greater span of asymmetries for given nuclear mass than for stable beams. 
Particularly important in the different areas, and simultaneously uncertain, is the density 
dependence of symmetry energy in uniform matter. As situation in the literature progresses, 
an increasing wider range of theoretical conclusions is getting made on that density depen- 
dence, as well as on some associated nuclear characteristics. In that situation, not unfamiliar 
from other contexts, rather than increasing claim statistics, it may be useful to take a step 
back and examine in detail the chain of implications following from the symmetry energy, 
to better understand the connections and to see how firmly any conclusions can be reached. 
We hope to make such a progress here within nuclear structure. 

Two main directions are advanced in our paper. One is in extending and combining 
symmetry-energy considerations from nuclear energy formula and from uniform nuclear mat- 
ter, within the nuclear energy functional of the Hohenberg-Kohn type {l| . Another direction 
is in carrying calculations of semi-infinite nuclear-matter within the Skyrme-Hartree-Fock 
approach \'}] . The framework developed in considering the energy functional is employed in 
analyzing the Hartree-Fock calculations. In the past, an analysis pertaining to the symmetry- 
energy, in the context of Hohenberg-Kohn functional, has been carried out by Farine |3|. 
Semi-infinite asymmetric matter has been examined before in the Hartree-Fock approach, 
in particular by Kohler {4, ^ and by Pearson et al. 6, [?, Q], and in the relativistic Hartree 
approach by Del Estal et al. [9]. Also, a number of Thomas- Fermi calculations for the 



half-infinite asymmetric matter have been done [1^, Q, Since the latest 

detailed Hartree-Fock studies by Pearson et al. [8'] have been done over two decades ago, 
there is potential for computational progress in these calculations, aside from our differing 
strategy in analyzing the results. 

In the following. Sec. [Ill reviews the basic nuclear energy formula, emphasizing the sym- 
metry energy. Section IIIII extends the symmetry-energy considerations to the context of 
a continuous nuclear Hohenberg-Kohn energy-functional. The limit of uniform matter is 
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discussed there too. In Sec. llVt details of our numerical calculations of half-infinite matter 
are presented and the results analyzed. Overall conclusions are presented in Sec. |Vl 



II. ENERGY FORMULA 
A. Elementary Formula 

The first basic context within which the nuclear symmetry energy is encountered is the 
empirical nuclear energy formula. In its most elementary textbook form, the formula con- 
tains just four macroscopic terms: 



Here, and Z are the neutron and proton numbers, respectively, A = + Z is the mass 
number and ay, a^, Oa and ac are the volume, surface, symmetry and Coulomb coefficients. 
The microscopic energy -Emic contains shell and pairing energy corrections. In spite of 
its simplicity, the energy formula is extremely powerful both in describing nuclear masses 
and in giving access the nature of nuclear interactions. When the microscopic corrections 
are disregarded, and the macroscopic coefficients are fitted to the data set on energies of 
over 3100 A > 10 nuclei, to yield ay = 15.3 MeV, as = 16.1 MeV, = 22.5 MeV and 
ac = 0.69 MeV, the rms deviation of the formula from data turns to be just 3.6 MeV, to 
be compared to the span of energies for the nuclei of over 2000 MeV. The nuclear energies 
are dominated by the volume term, responsible for nuclear binding, with binding energies 
reduced by the surface, symmetry and Coulomb terms. 

Regarding the nature of nuclear interactions, the dominance in E of the negative vol- 
ume term Ey, proportional to A, tells that the attractive nuclear interactions are short- 
range. The invariance of the symmetry term Ea with respect to neutron-proton interchange 
demonstrates that the nuclear interactions possess such an interchange symmetry, commonly 
termed charge symmetry. The positive value of suggests that neutron-proton interactions 
are more attractive than like-nucleon interactions; the Pauli principle effects contribute pos- 
itively to that coefficient as well. The surface term Es scales with A as in proportion to the 
surface of a nucleus with fixed volume per nucleon, represented as 47rrQ/3, or, equivalently, 



E{N, Z) = E^ac + E, 



mic 




av A + as A^'^ + a, 



(1) 
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with fixed density po = 3/(47rrQ). The reduction in the binding proportional to the surface 
may be thought of as one associated with the fact that nucleons at the surface experience less 
attraction than the nucleons in interior; the curvature of nuclear wavefunction at the surface 
also, in fact, contributes to that reduction. The form of the Coulomb term Ec is such as 
expected for a uniformly charged sphere, specifically (3/5) (Ze)"^ /Aneo R, where R = tqA^^^. 

Quantitatively, the fitted Coulomb parameter provides an estimate of the nuclear radius 
parameter, tq — (3/5) e^/(47reo etc) = 1-25 fm. With necessarily enhanced error, due to 
the raising of tq to cube power, the parameter further implies an estimate for the nuclear 
density po- Under changing nuclear density p, the nuclear energy should minimize at po, 
reaching the value of —ay per nucleon for N = Z and in absence of Coulomb interactions or 
boundaries. Given the surface term, one can estimate the nuclear surface tension, which is 
the energetic cost of creating the surface per unit area, or tension, a = dE/dH = as/47rrQ ~ 
0.8 MeV/fm^, again with a deteriorated expected accuracy due to the power of radius and 
due to combining of parameters. In the tension, E is nuclear surface area taking on the 
value of Atctq A?!"^ for a spherical shape. It is apparent that the energy formula can provide 
a wealth of information on nuclei. Naturally, this can be expected to extend to the changes in 
nuclear properties with changes in the relative neutron-proton asymmetry, r/ = (A^ — Z^j A. 

The success of the energy formula may be extended by modifying terms and incorporating 
new ones, with or without new parameters, including the microscopic terms. The terms may 
be guided by physics or ad hoc. The formula success can be perceived in at least two different 
ways. One is in terms of a better predictive power of the formula for net energy. In fact, 
by adding different terms and parameters, the rms deviation from fitted nuclear can be 



reduced by a factor of order of 5 [l7|, [l8|, [l9|, l20|, l2l|. Another measure of success is in 



enhanced access to fundamental nuclear properties. The advances in the two directions 
do not necessarily go hand in hand. Thus, obviously, adding an ad hoc term may throw 
off the value of a parameter that has a physical meaning, while helping though to reduce 
the rms error [2^. Even when a term of genuine physical meaning is added, problems can 
occur when many terms and parameters are present [22|. As an example, one can envision 
that the range of variation of A^ and Z for measured nuclei may not be suffice to make 
two terms in an energy formula distinct, with different parameter combinations yielding 
comparable agreement with data and a parameter expected a •priori to be obscure, eating into 
a major parameter. Another example might be a term introduced to improve the description 
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of one N-Z region (e.g. low A) that gets its parameters determined, in an unguided fit, 
by the prevalence of nuclei in another N-Z region (such as high A). The problems with 
multiparametric fits are known from other areas and the developed remedy is to introduce 
as many prejudices [23] into the fit as possible. The prejudices may, in fact, include the 
qualitative features of fitted functions, following from physical considerations, common-sense 
limits on region for parameter determination and the use of constraints on the parameters 
from auxiliary investigations. The latter constraints, in the case of an en ergy formula, could 



e.g. involve information on the proton and neutron densities, and p„ 20|, 124, l25 |. 



B. Surface Symmetry Energy 

Examining the symmetry term in the basic formula ([T]), one can notice that this term 
has a volume character, i.e. it changes in proportion to A, like Ey^ when N and Z are 
changed by same factor. One question which arises, when addressing the symmetry energy, 
is whether there is a conceptual need for adding a surface symmetry term to the energy 
formula. Another question is how such term should enter that formula. In the literature, 
there are at least two ways of introducing such a term, with a surface symmetry parameter 
having different meanings in the two formulations. 

Regarding the first question, the volume energy becomes less negative, when magnitude 
of nuclear asymmetry increases from zero, at a fixed A. With this, less work should be 
required for developing a surface, since there is less lost binding then to compensate for 
compared to zero asymmetry. Thus the surface tension should drop with increase in asym- 
metry magnitude. The quantity associated with asymmetry, shared by systems in contact in 
macroscopic equilibrium, such as the surface region and interior, is the asymmetric chemical 
potential 

^'^ = d{N - Z) = 2^^- - ■ 
Under charge symmetry of nuclear interactions, with Coulomb ignored for the moment, the 

nuclear tension should depend only on square of pa and, following the above, behave for 

small asymmetries as 

dEs 2 /Q^ 

^ = ^ = ^0 -Z^/^a, (3) 

where (Tq = as/4:iTrQ is the tension for symmetric nuclear matter and u is some positive 
constant. When the tension depends on asymmetry, so must the surface energy Es, cf. ([3]). 




FIG. 1: Surface asymmetry can be understood in terms of different radii for proton and neutron 
distributions. 

The next question is of the form of an energy formula which incorporates that dependence 
and, of course, of the consequences of that dependence. 

The inverse Legendre transformation for asymmetry, from fia to N — Z , is given by 

9$ 



N-Z 



(4) 



(5) 



where 

$ = iia{N -Z)-E = fia{N -Z)-Ev-Es. 
Upon inserting of ([5]) into (jlj), we find 

d{N-Z) _^dE _,(dEy ^ dEs\ , 

In the equation above, N — Z, E, Ey and Es are all extensive, i.e. they grow, in a character- 
istic fashion, as the system size grows. The equation implies that, when the surface energy 
depends on fia, the net asymmetry partitions into volume and surface portions, Ny — Zy 
and Ns — Zs, satisfying respectively 



d{N-Z)y^S_^^dEy,S 



(7) 



(9/ia d/j^l 

The surface asymmetry may be transcribed onto a difference of radii for neutron and proton 
distributions, which is the basis of the droplet model (to be discussed) of nuclei and is 
schematically indicated in Fig. [H 

Given that, at low asymmetries, due to charge symmetry, the energies must depend 
quadratically on asymmetry, as does the tension in Eq. Q, the second derivatives of energy 
in Eqs. and ([7]) can only depend on A. In consequence, at low asymmetries, also the 
first derivatives of — Z in those equations must also depend solely on A. Accounting for 
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the dimensions and for scalings with A, we can write 

A v42/3 
Nv - Zv = — ^7 yWa and Ns - Zs = -—g fia , (8) 

where and af are constants characterizing the nuclear volume and surface, respectively, 
with the dimension of energy. Sum of the two asymmetries gives us the relation between Ha 
and net asymmetry: 

N-Z = Nv-Zv + Ns-Zs = ^(4 + ^) ■ (9) 
From Eqs. and (IHl), we further find 

Ey = E'y + ^f^l = -ayA + al ~ , (10) 

and 



Es - + —^fi^-asA ' — , (11) 



where -Ey^ are the volume and surface energies for symmetric N = Z matter. It should be 
mentioned that there is no error in the sign difference in symmetry terms between Eqs. (ITT]) 
and Whereas the volume and surface energies increase with asymmetry, the surface 
tension decreases. The derivative of surface energy with respect to surface area, for the 
surface tension a, must be taken at constant Ns — Zs, which yields i>~^ = 16 



2 „s 



On adding up the energy contributions ffTOj) and ffTTl) . we get for the net energy 

.^.»4(A,^)._...,....3,M, 

^ \"a "a / -y + 

where we have eliminated fia in favor of N — Z using (jl]). The nuclear energies above, 
quadratic in asymmetry, exhibit analogy to the capacitor energy in terms of electric charge. 
While asymmetry is analogous to the electric charge within this analogy, the chemical poten- 
tial is analogous to the electric potential. The coefficients of proportionality between volume 
and surface asymmetries and chemical potential in ([8]) are analogs of capacitance, with the 
two capacitances proportional, respectively, to the volume and surface. For connected ca- 
pacitors, the potentials are equal. The net symmetry partitions itself in (Q, between the 
volume and surface, in proportion to the capacitances. Finally, net energy for connected 
capacitors can be represented in f|T2l) either in terms of net potential squared multiplied by 
net capacitance or in terms of net asymmetry squared divided by the net capacitance ^. 



^ Beyond the quadratic approximation for the energy in asymmetry, the capacitance may be defined 26 1 in 
terms of the derivative of asymmetry with respect to the chemical potential, i.e. the l.h.s. of 
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Upon adding the Coulomb and microscopic contributions, we now arrive at the energy 
formula with a mass- dependent symmetry coefficient 

E{N, Z) = -ayA + as A'/' + ^ (iV - Z)^ + ^ + E^,, , (13) 



where 



In the limit of large A, the asymmetry gets primarily stored within the volume and then the 
coefficient aa{A) tends towards a^. On the other hand, in the limit of small A, the storage 
of asymmetry gets shifted to the surface and then the ratio aa{A)/A approaches a^/A^^^. 

Although the surface symmetry effects can be important in learning on symmetry energy, 
as far as the net nuclear energy is concerned, those effects just correct the symmetry energy 
which itself is a correction to the leading term —ay A term in an nuclear energy. Additionally, 
since most of the measured asymmetric nuclei are heavy, those nuclei can provide little 
variation for A^/"^ in the symmetry-energy term. In consequence, the basic energy formula ([1]) 
can provide a ffi to the measured nuclear energies that is quite satisfactory without any 
surface symmetry modification. 

The coefficient af of surface symmetry energy is generally of interest as a fundamental 
quantity and because, in the context of the changing density at nuclear surface, it can provide 
access to the density dependence of symmetry energy in bulk nuclear matter. However, 
obviously, since af appears within a correction to a correction in the energy formula, it 
can be difficult to learn about that coefficient by ffiting nuclear energy data with a formula 
where further types of secondary corrections may need to be included. The latter corrections 



may compete against each other and against the surface symmetry energy 22|, |25| . This can 
be especially true for blind ffis. The potential secondary corrections include those to the 
asymmetry- dependeni Coulomb term and those to the a,sym.m.eiiy- independent part of the 
formula. The Coulomb energy e.g. also affects the surface asymmetry. Corrections to the 
asymmetry-independent part of the formula are important, possibly against expectations, 
because the grouping of nuclei around the hne of stability, in the (iV, Z)-plane, introduces 
mass- asymmetry correlations for the ffited data. 

In this and in a subsequent paper, we shall try to understand the interplay between the 
dependence of bulk symmetry energy on density, density distributions in nuclei and the 
symmetry coefficient aa{A). On one hand, we shall try understand a connection between af 



and the density dependence of symmetry energy. On the other hand, we shall try to deter- 
mine aa{A) from nuclear data, refraining from a blind fit and without insisting on validity of 
Eq. (IT^ . Finally, we shall seek additional information on symmetry energy within density 
distributions. 



C. Droplet Model 

A result of the forr a f|l^ . for the symmetry coefficient, has first appeared in the literature 
in the droplet model 27| by Myers and Swiatecki. The droplet model relies on the assump- 
tion that the neutron and proton surfaces shift relative to each other as nuclear asymmetry 
changes. The surface energy per the elementary area of An Tq depends in the model on that 
shift, Dnp, and on the relative asymmetry in nuclear interior, riv, according to: 

as + H [ — +2Pr]v — --Gr]y. (15) 



In the above, H, P and G are model constants and we have dropped a higher-order curvature 



term employed in [27|. (Strictly, the surface energy density Es/Tj is identified in j27| as 
the surface tension a, but these two are different when they depend on asymmetr y, c f. 
Subsection IIIBj Eq. (fT5|) is consistent with the use of the formula by the authors of |27|.) 



When a higher-order Coulomb correction included in [271] is disregarded, the symmetry 
coefficient emerges such as given by Eq. ([HI), with the surface symmetry coefficient given 
by 

5 4^ 8 H^G 

where Q is an auxiliary constant. 



Following 2J] and Subsection III Bl Eq. (fT^ must emerge in the macroscopic limit at 
large A, for an energy analytic in asymmetry, provided underlying model assumptions rep- 
resent short-range interactions obeying charge symmetry. Later in the paper, we shall en- 
counter a greater richness of changes with asymmetry in the nuclear surface region for dif- 
ferent Skyrme interactions, in the Hartree-Fock calculations, than envisioned in the droplet 
model. However, all those calculations will be consistent with the symmetry energy being 
quadratic in asymmetry in the macroscopic limit, with the symmetry coefficient conforming 



with (ITI) at large A. 
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The droplet-model energy formula eventually acquired [l7| as many as 38 independent 
parameters that either got set by fitting nuclear energies or using auxiliary information. 
The general difficulty in settling on the value of af in an energy formula is illustrated 
by the fact that this coefficient has risen by over a factor of 2 within the history of the 
droplet model, from af = 7.1 MeV [27| to 15.7 MeV [28|. At the same time, the volume 
symmetry constant of the model has varied within the range of only about 25%, from 



al = 28.1 MeV 



27| to = 36.5 MeV 29|]. A still wider range of surface symmetry 



coefficients is found when considering other formulas in the literature. Thus, e.g. the fits 
by Pomorski and Dudek 19|], who have emphasized asymmetry effects, give rise to the 
coefficient values of = (fi;voi^voi)^/('tsurf ^surf 

) = 21.3 MeV and 50.1 MeV, in terms of these 
authors' notation, depending on the variant of the drop formula they employ in describing 
nuclear ground-state energies and fission barriers. Depending on the formula variant, the 
fits, at the same time, produce the values for the volume symmetry coefficient of either 



= -K^oi b^^i = 28^MeV or 25.4 MeV. Notably, the droplet formula [17[ and the Lubhn- 
Strasbourg formula [19] produce comparable rms deviations from the measured energies and 
barriers. Similarly wide variations of the surface symmetry coefficient, of nearly an order of 
magnitude, will be found in Sec. IIVI for the effective Skyrme interactions employed in the 
literature. The volume symmetry coefficients will, again, exhibit relatively less variation. 
The values of af will be found to be correlated with the density dependence of symmetry 
energy. 



III. ENERGY FUNCTIONAL 
A. Hohenberg-Kohn Functional 

Here, we shall consider the nuclear energy functional in terms of neutron Pn{f) and 
proton Pp{r) densities. Given a system hamiltonian in terms of nuclear and Coulomb com- 
ponents, 

H = -ffnucl + He , (17) 

different densities Pn{f) and proton Pp{r) and different values of the functional E{pn,Pp) will 
result, when minimizing the expectation value of the hamiltonian in the presence of different 
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external potentials [ll, T4(r) and Vp{r): 

{H + Vn + Vp) = E{pn, pp) + Jdr (K(r) p„(r) + Vp{r) pp(r)) . (18) 

The functional E, defined in this way, is specified for integer nucleon numbers, N and Z, 
only. In finding the ground-state densities and energy, the functional E is minimized with 
respect to variations in p„ and Pp. In the context of an energy formula, we will actually 
be interested in the functional E averaged out and smoothed out across discrete and Z, 
defined then also for fractional nucleon numbers and made analytic in the nucleon densities. 

The smoothed-out functional E will be partitioned into the parts associated with nuclear 
and with Coulomb interactions. The nuclear part will be further partitioned into the energy 
of symmetric matter and a correction associated with asymmetry: 

Eipn, Pp) = ^nucl(P, Pnp) + ^c(P, Pnp) = ^o(p) + ^a(P, Pnp) + ^c(P, Pnp) , (19) 

where p = p„ + p^, p„p = p„ - Pp and Eo{p) = ^nuci(p, 0). 

At Pnp = 0, the symmetry energy Ea vanishes by definition. Given that energies of nuclei, 
corrected for microscopic and Coulomb contributions, minimize at a lowest asymmetry, the 
energy Ea must be quadratic in asymmetry at low asymmetries: 

Ea{p,pnp) = j dridr2pnpiri)S{p,ri,r2) pnp{r2) + 0{r]^) . (20) 

Because of the minimum, S must be a positive-definite symmetric bilinear integral operator, 
with regard to its arguments Vi and This operator represents an analog of the potential 
matrix in electrostatics, dependent on conductor described in terms of p. For the purposes of 
assessing later the accuracy of our considerations, we have indicated the principal presence 
of the subsequent terms of fourth order in asymmetry and higher. At other times, we may 
omit underscoring the presence or order of any subsequent terms in an expansion, beyond 
the directly considered order. 

Minimal energy E at given N and Z will be achieved for densities p and pnp which satisfy 
the equations 

. (21) 



6p{r) 

and 

6E 
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where jj, and /ia are Lagrange multipliers associated with the constraints 

jdrpir)=A, (23) 

and 

J drpnp{r)=N-Z, (24) 

respectively. The Lagrange multipliers coincide with the corresponding chemical potentials, 
which are derivatives of energy E that gets minimized, fi = BE /OA and /ia = dE/d{N — Z), 
hence the notation for multipliers. Thus, e.g. upon differentiating the energy with respect 
to mass number, 

d-E f.SE dpir) f dfir) 

we find that the l.h.s. coincides with fi introduced as a Lagrange multiplier in ( 12T1) . Since 
the nucleonic chemical potentials are equal to 

Atn,p = ± /ia , (26) 

the system can spontaneously emit nucleons, either neutrons or protons, if 

\^^a\>-^i■ (27) 

If we consider the symmetric term Eq alone, we can expand it around the density po{r) 
that minimizes Eq amongst the densities that meet the condition fl23l) at a prescribed A: 



Eo{p) = Eoipo) + fio J drAp{r)+ J dn dra Ap(ri) /C(ri,r2) Ap(r2) + C ((Ap)^) 
= EoiA) + J dridr2Ap(ri)/C(ri,r2)Ap(r2) + 0((Ap)3) , 



(28) 



where Ap{r) = p{r) — poir), po = dEo/dA, /C is a positive-definite operator and / dr Ap = 
holds amongst the densities that meet the condition (l23l) . 

Note that, depending on context, we may use the symbol po to denote either the function 
of position that minimizes the energy of a finite or semi-finite symmetric system, in the ab- 
sence of Coulomb interactions, or to denote single density value that minimizes the energy of 
symmetric uniform matter. Use of a position argument for the quantity will obviously imply 
the function. The distinction may be, otherwise, made by directly naming the quantity. 
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B. Densities and Generalized Symmetry Coefficient 



In the following, we shall suppress the Coulomb term Ec and we shall explore what con- 
sequences the discussed features of the functional E may have on the asymmetry dependence 
of nuclear densities and of nuclear energy. 

Thus, at given N and Z, for low asymmetry \N — Z\, the minimal energy may be generally 
represented as 

E{N, Z) = Eo{A) + ^ (iV - Zf + 0{r^') , (29) 

where aa{A) is a general A-dependent symmetry coefficient. From fl29|) . the asymmetry 
chemical potential, in terms of the coefficient, is 

= iN-Z) + Oiv') , (30) 

while the symmetric chemical potential from fl29|) is 

^, = ^,o + {N-Zr^^ + Oiv'). (31) 

When the Coulomb term is suppressed, Eqs. (12T|) and (l22l) for the densities produce 

fia = 2j dri Sip, r, ri) p„p(ri) + 0{r]^) , (32) 

and 

= 1^ + / dri dr, p„,(rO ^^^^^^^'^j''^^ Pn^) + Oiv') ■ (33) 
An operator inverse to S, 

j dr2S{p,ri,r2)S'\p,r2,r3) = J dr2S'\p,ri,r2) S{p,r2,r-i) = 6{ri - r^^) , (34) 

should exist as S must have positive eigenvalues only. Upon applying the inverse operator 
to both sides of (|32|) . we can formally solve that equation to get 

Pnp(r) = dr,S-\p,r,r,) + 0{r^') . (35) 

Upon integrating both sides of (135!) over space, we find for the asymmetry chemical potential 

J dridr2S-\p,ri,r2) 



As a consequence, the nucleon density difference from fl35|) can be represented as 

p.,(r) = (iV - Z) /dn^^^(p,r,r,) ^ ^ ^ 
J dridr2 5-i(p,ri,r2) 
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Upon substituting (1371) into the symmetry energy (l20l) . we get 



Ea = y 7 + 0{r^') . (38) 



Substitutions into the equation fl33|) for the net density, next yield 

''» = iT\ + (^-^''7TT(77rWv (39) 

5p(r) 5p(r) Vj dridr2 5-i(p,ri,r2) A J 

where we have used 6A/S p(r) = 1 and the identity 

y dr2 ^5(p,n,r2) ^^^^^ + ^^^^^ 5 (p,r2,r3)J=0, (40) 
following from differentiating Eq. ( IMl) side-by-side . With po(5') representing the solution of 

5p(r) 

we find from ( |39l) that the net density for a finite asymmetry differs from that for symmetric 
matter only up to second order in asymmetry, 

p(r) = po(r) + 0{v') , (42) 



j drdriS-\po,r,ri). (43) 



with the second order brought in by the second term on the r.h.s. of ( l39l) . Upon comparing 
the expanded energy in ( l29l) to the combination of energies Eq from ( l28l) and -E^ from ( l38l) . 
with the density po{r) from fH2l) . we arrive at a result for the generalized symmetry coefficient 
in terms of the operator S, 

A 

aaiA) 

Given the weak dependence on asymmetry for the net density, stated in Eq. fH2|) . we 
conclude, from (!35|) or (1371) . a similarly weak dependence on asymmetry for the density 
difference pnp scaled by pa or {N — Z) (or another quantity odd in asymmetry). In that 
context, we introduce an isovector density changing weakly with asymmetry, that will be 
normalized either locally or globally. Normalized locally, in terms of the intense chemical 
potential, the isovector density is defined as 

Pair) = ^ p„,(r) = ^ M!l + ^^^2^ (44) 

Pa aa{A) 1] 

where the r.h.s. follows from (1301) . For an extended system at low asymmetry, the dimension- 
less ratio normalizing pnp on the r.h.s. of (j44l) can be identified with volume asymmetry r]v, 
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cf. Eq. dH]). The normalized chemical potential ^a/'^o)^ differs from that asymmetry only by 
higher order terms and may be termed effective volume asymmetry: 



^'v^V^ = ^,^ + ^iv') =Vv + 0{v') . (45) 
2a^ aa{A) 



In the limit of symmetric matter, the isovector density becomes 



dpnp{r) 



dpnp{r) 



(46) 

ri=0 



Normalized globally, to integrate to a mass number A, the isovector density pa is defined as 
PA{r) = — p„,(r) = ^ p,(r) = ^ p,(r) + 0{v') . (47) 



In the limit of symmetric matter, the density normalized that way becomes 

P» = ^ • (48) 

From fl35|) , it follows that the isovector density pa may be expressed in terms of the operator S 

as 

Pa{r) = al j <\T,S-\p,r,r,) + 0{ri^)=pl{r) + 0{ri^), (49) 
with the 7] = limit of 

Plir) = al j dnS-\p\r,n). (50) 

From fH3l) . it follows that the symmetry coefficient is related to the locally normalized density 
with 

The result ( 1491) for the isovector density is a counterpart of the result (142|) for the isoscalar 
net density. The obvious important implication of those results is that the proton and 
neutron densities should be describable, up to second-order in asymmetry and up to Coulomb 
corrections discussed in the next paper, in terms of two densities universal across an isobar 
chain, 

Pn,v{r) = \ {p{r) ± p^(r)) = 1 (p(r) ± p,(r)) = 1 (po(r) ± ^ p°(r)) + 0{r,^) . (52) 

As is apparent, the isovector density is tied to the symmetry coefficient. We shall investigate 
quantitatively the relation between the densities p and pa for different interactions in half- 
infinite matter in the next section. It should be mentioned that, while pa represents an 
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isovector form factor, it is, at this point, an isoscalar with regard to its transformation 
properties in isobar space. Using analogy with spin, that density might further be termed 
an isospin susceptibihty. 

One issue that Eq. flS21) raises is that of the importance of higher-order terms in asymmetry 
expansion. That importance needs to be considered not just globally but locally as well. 
In the next section, we shall see that pa is typically larger than p outside of the nuclear 
surface. Since nucleon densities must be positive, the higher-order terms must become 
significant, at a given r^y, at distances outside of the surface where 

lw|Pa(r)>Po(r). (53) 



C. Uniform Matter 



In uniform matter, the energy is an integral of constant energy per unit volume e, 

/E E 
dre{pn,pp) = Ve, with e = — = p — , (54) 



and 



| = f (p) + f (p,^)-f (p) + 5(p)r/^ (55) 



where S" is a density-dependent symmetry coefficient. That last coefficient is on its own 
commonly termed the symmetry energy. An obvious test of the usefulness of the quadratic 
expansion of symmetry energy is whether the approximate equality, 

§(P,1)^5(P), (56) 

holds at different densities, particularly at p = po- 

To the extend that the approximation of the r.h.s. of fl55l) holds js^, the functions -^(p) 
and S{p) allow for the determination, at any asymmetry in uniform matter, of nuclear energy 
and nuclear pressure, 

(E\ 2d (E,\ ^ 2 2d5 

and nucleonic chemical potentials, 

P.„ = P^^^J+^ = P±P., (58) 
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where 

d /E\ E P E 

and 

d fE 



The pressure and chemical potentials are important in calculating neutron-star properties 
and in simulating supernova explosions 3ll ]. 

Around pQ, it is common to expand both Eq/A, 



_(p) = -a^ + (p-po) +..., (61) 

and 5, 



= «r + (P - Po) + (P - Po)^ + . • • ■ (62) 

Notation for the constant terms in the expansion is consistent with that in the energy 
formula, Sec.[Tll Upon employing the expansions fl6Tl) and (|62|) . the energy per nucleon (|55|) 
may be rewritten, around its minimum, in the form 

— = -av + a^7] + ^ [p - Pmin) + ■ ■ • , (63) 
^ -LoPmin 

where we retain an accuracy of the order of 0{ri'^) in the coefficients. Within that accuracy 
the density which minimizes the energy is 

3L 

Pmin(^) - Po - Po ^^ (64) 

and the effective incompressibility is 

i^eff(r7) = K+ (i^sym " 6L) 7/^ . (65) 

In self-bound matter, on one hand, the pressure vanishes, P = 0, and, on the other, the 
nucleonic potentials are negative. When assuming a density close to normal, from ( I58l) -( 1MI) . 
one finds that the magnitude of bulk drip asymmetry rjy, above which nuclear matter emits 
nucleons, is approximately 

In gravitationally bound neutron-star matter, due to the fact that the energy of symmetric 
matter minimizes at po, the constant L from the symmetry-energy expansion (162|) . largely 
determines the nuclear contribution to pressure in neutron-star matter around po, 

P-v^p'j-p-^p^'n'- (67) 
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Table U shows a selection of bulk nuclear parameters for the majority of Skyrme interac- 
tions employed in the literature in nuclear mean-field calculations. The force constants for 



most of those interactions have been compiled by Jirina Stone |32l. |33[|. Basic information on 
the determination of Eq/ A and S of uniform matter for those interactions will be provided 
later in this section; more thorough explorations of the properties of uniform matter may 



be found elsewhere, e.g. in 



32l |. On the other hand, procedures for the determination of 



surface parameters and values of those parameters for the interactions will be discussed in 
the following section. 

Variations in the bulk parameter values in Table [T] can serve as illustration of the status 
of the specific parameters. Thus, there is a general consensus in the literature regarding 
ay and po? and the respective values in the table, across the listed interactions, are ay = 
15.9 ± 0.3 MeV and po = 0.158 ± 0.005 fm"^. Regarding variation of the energy of symmetric 
matter with density, recent analyses of the energies of giant isoscalar monopole resonances 
point 3J] to nuclear incompressibilities within the range of = (230 — 250) MeV, and fewer 
than 30% of all the interactions in Table [T] give rise to incompressibilities definitely outside 
of that range. Some of the interactions characterized by high i^"-values outside of the above 
range have been, incidentally, purposely constructed to explore the effects of high i^'-values 
on properties of finite nuclei. 

Regarding the dependence oi E/A on rj, the values of expansion coefficients for normal 
Po-matter in Table [H a^, are lower by just ~ 1 MeV than the values of symmetry energy 
per nucleon at ?7 = 1, Ea/A{po, 1). This indicates that the quadratic expansion of energy in 
asymmetry is fairly accurate near po for the Skyrme interactions. The next quartic coeffi- 
cients in the expansion oi E/A with respect to asymmetry must be positive for the Skyrme 
interactions and be of the order of 1 MeV. The inferred magnitude of quartic terms indi- 
cates that a term breaking charge invariance, that might be potentially introduced into E/A, 
could be more important around po than the quartic term, for typical rj of interest. 

The relative variation across Table [T] is far greater for = S{po) than for ay. The low- 
est values, of the order of 23 MeV, are close to the Oa-value inferred from the basic 
energy formula ([1]) from Sec. [Tll while some of the high values in the table actually ex- 
ceed 37 MeV. An even greater degree of variation is found in the table for the constants L 
and -ft'sym; characterizing the p-dependence of S. Thus, the derivative-constant L takes on 
positive and negative values in the table, within principally as much as 3 orders of magni- 
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FIG. 2: Symmetry energy 5 as a function of density for sample Skyrme interactions characterized 
by widely varying derivative-parameters L. 

tude; the average value of L there is 35 ± 33 MeV. The situation is found to be similar for 
i^^sym that takes on both negative and positive values within 2 orders of magnitude. Func- 
tions S{p) for the Skyrme interactions representing two more typical and two more extreme 
L- values are shown in Fig. [2l 

While the parameters L and ii'sym vary widely between interactions, they do so in a fairly 
correlated manner. The relatively tight correlation is illustrated in Fig. [3] displaying results 
for the Skyrme interactions, in the plane of those parameters scaled by a^. The correlation 
may be understood as due to the fact that just one physical scale of po governs the low- 
density variation of nuclear energy, limiting possible independent variations in the first and 
second derivatives of S at po- If we were to approximate 5* (p) -functions for the Skyrme 
interactions, such as those shown in Fig. [21 with parabolas, S{p) oc p(l + cp), i.e. extending 
the expansion fl62|) to the whole subnormal region of densities, we would arrive at the 
following relation between the scaled parameters: 

^«6-^_i8. (68) 

That relation, represented by the solid line in Fig. [3], qualitatively explains the correlation for 
the interactions in the figure. However, a parabola produces an incorrect linear behavior of S 
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FIG. 3: Correlation between the symmetry-energy parameter-ratios Kf^y^/a^ and L/a^ . The cir- 
cles represent Skyrme interactions, while the lines represent the modified power parameterization 
S{p)^p'^{l + cp). 

at low p, where 5* is expected to be dominated by the Fermi-energy contribution proportional 
to p^^^. An attempt to improve the approximation to S, using S{p) oc p'^(l -|- cp), leads to 
the relation 

T<r T 

(69) 



represented by the dashed line in Fig. [3l for z/ = 2/3. It is apparent that the v = \ and 
z/ = 2/3 relations practically bracket the correlation between the L and -ft'sym parameters, 
in spite of the fact that S for the Skyrme interactions differs from either of the forms 
used in the approximations. Practical implication of the observed correlation is that any 
constraining of either of the two parameter will produce constraints on the likely values of 
the other parameter. 

With ay ~ 16 MeV, and with the average volume symmetry parameter for the Skyrme 
interactions being ~ 30 MeV, the magnitude of drip asymmetry expected for the Skyrme 
interactions from the simple expression (!66l) is ~ 0.27, which turns out to be close to the 
typical magnitude of r/y in Table [B The simple formula implies further an anticorrelation 
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between the bulk asymmetry coefficient and the drip asymmetry. Indeed, the lowest and the 
highest drip asymmetry values in Table [T] are associated with some of the highest and the 
lowest, respectively, volume coefficient values in the table. Overall, the magnitudes of 
rjy in the table may seem to be low. However, these values pertain to nuclear interiors. For 
a given Skyrme interaction, due to asymmetry effects in the surface, nuclei may generally 
have a net asymmetry t] larger in magnitude than rjy and still be stable. 
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TABLE I: Bulk nuclear properties for different Skyrme interactions with force constants given in the indicated references. 



Pq ay m* /m K L 
(fm-3) (MeV) (MeV) (MeV) (MeV) 



K 

(MeV) 



f (Po,l) 
(MeV) 



as 
(MeV) 



(MeV) 



do AR^ 
(fm) (fm) 



Ae-R (i° Ref. 
(fm) (fm) 



0.1553 -15.99 0.911 370.3 29.24 1.2 -461.8 30.27 0.325 17.36 26.00 0.438 0.432 

0.1482 -15.96 0.580 340.8 34.14 50.1 -265.1 35.75 0.274 19.38 17.54 0.517 0.718 

0.1453 -15.85 0.763 355.3 28.16 9.9 -393.7 29.38 0.337 18.54 21.77 0.486 0.518 

0.1507 -16.57 0.789 371.8 32.66 29.3 -367.7 33.60 0.299 20.06 17.70 0.492 0.653 

0.1509 -15.96 0.471 324.5 31.22 63.5 -136.7 33.22 0.312 18.87 13.39 0.529 0.779 

0.1551 -16.05 0.383 305.7 32.83 96.1 24.2 35.32 0.312 19.16 10.27 0.553 0.941 

0.1435 -15.75 0.949 363.6 26.88 -7.3 -471.3 27.86 0.354 18.10 26.27 0.466 0.428 

0.1434 -15.79 1.001 366.4 26.96 -10.1 -488.9 27.89 0.354 18.11 27.43 0.462 0.415 

0.1476 -15.40 0.602 333.3 24.89 28.2 -236.7 26.25 0.385 14.21 17.58 0.431 0.477 

0.1610 -15.98 1.000 236.2 32.02 56.2 -134.9 32.73 0.305 18.26 14.60 0.556 0.799 

0.1610 -15.94 1.000 235.8 32.00 56.2 -134.7 32.71 0.305 18.01 14.71 0.552 0.794 

0.1610 -15.94 1.000 235.8 31.50 55.3 -132.1 32.21 0.311 17.78 15.33 0.547 0.776 

0.1590 -15.95 1.000 235.5 35.46 94.1 -24.5 36.16 0.284 18.17 11.57 0.558 0.986 

0.1640 -16.00 1.000 201.7 37.01 98.5 -25.0 37.72 0.274 18.13 10.91 0.604 1.084 

0.1609 -15.96 1.000 236.0 29.97 30.9 -211.6 30.67 0.321 18.22 18.10 0.555 0.658 

0.1606 -15.94 0.833 235.7 29.52 31.1 -209.9 30.72 0.326 18.12 18.21 0.559 0.650 
Table Continued on Next Page. . . 



0.433 0.461 [35] 

0.761 0.439 [35] 

0.509 0.474 [36] 

0.717 0.412 \37] 

0.906 0.404 [36] 

1.231 0.368 [3^ 

0.404 0.496 [36] 

0.388 0.500 [37] 

0.554 0.393 \38\ 

0.834 0.450 [39] 

0.828 0.450 [39] 

0.782 0.468 [39] 

1.171 0.390 [39j 

1.283 0.400 [39] 

0.630 0.508 \39} 

0.617 0.516 [39] 
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Name 



Po 



ay m*/m K 



sym 

(fm-3) (MeV) (MeV) (MeV) (MeV) (MeV) 



(MeV) 



4 



as 
(MeV) 



(MeV) 



(fm) (fm) 



Ae-R (i° Ref. 
(fm) (fm) 



SkT8 

SkT9 

SkTK 

SkM 

SkMl 

SkMP 

SkMs 

SKa 

SKb 

SGI 

SGII 

RATP 

SkP 

E 

Es 

Gs 

Rs 



0.1607 -15.94 0.833 235.7 29.92 33.7 -187.6 30.77 0.322 18.13 17.88 0.559 0.657 0.637 0. 

0.1603 -15.88 0.833 234.9 29.76 33.8 -185.7 30.60 0.323 17.79 17.81 0.555 0.654 0.636 0. 

0.1681 -16.70 0.611 253.3 35.57 41.6 -221.9 37.35 0.275 19.78 19.46 0.579 0.696 0.685 0, 

0.1603 -15.77 0.789 216.7 30.75 49.4 -148.9 32.07 0.314 16.85 14.84 0.558 0.756 0.789 0, 

0.1603 -15.77 0.789 216.7 25.17 -35.3 -389.0 26.53 0.414 17.46 59.67 0.575 0.180 0.161 0, 

0.1570 -15.56 0.654 230.9 29.89 70.3 -49.8 31.26 0.336 16.58 12.06 0.549 0.850 0.950 0. 

0.1603 -15.77 0.789 216.7 30.03 45.8 -156.0 31.39 0.322 17.46 14.48 0.575 0.765 0.790 0. 

0.1554 -15.99 0.608 263.2 32.91 74.6 -78.5 34.56 0.301 18.77 12.51 0.561 0.877 1.013 0. 

0.1554 -15.99 0.608 263.2 23.88 47.6 -78.5 25.53 0.460 18.77 10.27 0.561 0.768 0.895 0. 

0.1544 -15.89 0.608 261.8 28.33 63.9 -52.0 29.62 0.365 17.48 12.76 0.531 0.770 0.856 0, 

0.1583 -15.59 0.786 214.7 26.83 37.6 -146.0 28.10 0.365 16.10 15.11 0.538 0.660 0.679 0, 

0.1598 -16.05 0.667 239.6 29.26 32.4 -191.3 30.80 0.332 18.68 17.95 0.579 0.658 0.621 0, 

0.1625 -15.95 1.000 201.0 30.00 19.7 -266.7 31.33 0.316 18.18 18.02 0.607 0.670 0.631 0, 

0.1591 -16.13 0.868 333.4 27.66 -31.2 -570.8 28.83 0.359 18.94 37.62 0.485 0.321 0.280 0. 

0.1628 -16.02 0.839 248.6 26.44 -36.8 -457.9 27.72 0.390 18.21 52.37 0.543 0.222 0.191 0. 

0.1576 -15.59 0.784 237.3 31.37 94.0 14.0 32.61 0.336 16.07 10.10 0.512 0.929 1.191 0. 

0.1577 -15.59 0.783 237.4 30.58 85.7 -9.1 31.82 0.340 16.06 10.59 0.512 0.888 1.106 0, 
Table Continued on Next Page. . . 
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508 [39 



512 
460 [41] 



729 [42 
423 [43 



470 [44] 
409 \5} 
408 [5J 
424 [45] 



464 [45 
534 [46 



533 [43 

576 [48] 



689 [48 
357 [48 
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TABLE I - Continued 



Name 



Po 



ay m*/m K 



sym 

(fm-3) (MeV) (MeV) (MeV) (MeV) (MeV) 



(MeV) 



(MeV) 



(MeV) 



(fm) (fm) 



Ae-R (i° Ref. 
(fm) (fm) 



T 
Z 
Zs 

Zss 

SkSCl 

SkSC2 

SkSC3 

SkSC4 

SkSC4o 

SkSC5 

SkSC6 

SkSClO 

SkSCll 

SkSC14 

SkSC15 

Skyrmelp 

MSkA 



0.1613 -15.93 1.000 235.7 28.35 27.2 -206.8 29.06 0.343 17.74 22.64 0.547 0.587 0.476 

0.1589 -15.97 0.842 330.3 26.82 -49.7 -657.9 27.97 0.391 17.74 51.51 0.465 0.213 0.199 

0.1630 -15.88 0.783 233.4 26.69 -29.3 -401.6 28.01 0.371 17.00 46.62 0.534 0.233 0.217 

0.1625 -15.96 0.775 234.9 28.80 -4.5 -332.7 30.14 0.330 17.32 29.34 0.541 0.406 0.372 

0.1607 -15.85 1.000 234.6 28.10 0.2 -312.1 28.81 0.338 17.28 26.64 0.536 0.448 0.401 

0.1607 -15.90 1.000 235.2 24.74 11.0 -228.3 25.44 0.401 17.45 19.73 0.540 0.515 0.477 

0.1606 -15.85 1.000 234.5 27.01 0.8 -296.3 27.71 0.355 16.88 27.39 0.529 0.434 0.375 

0.1606 -15.86 1.000 234.8 28.80 -2.1 -329.6 29.51 0.328 17.37 28.15 0.538 0.439 0.389 

0.1606 -15.84 1.000 234.5 26.98 -9.6 -337.8 27.69 0.357 17.35 30.46 0.538 0.385 0.337 

0.1606 -15.85 1.000 234.5 30.99 -6.9 -375.2 31.69 0.300 17.28 32.03 0.536 0.415 0.368 

0.1607 -15.92 1.000 235.5 24.57 11.0 -226.3 25.28 0.406 17.59 19.56 0.543 0.518 0.478 

0.1607 -15.96 1.000 235.9 22.83 19.1 -172.8 23.53 0.457 17.75 15.61 0.546 0.577 0.556 

0.1606 -15.86 1.000 234.8 28.80 -2.1 -329.6 29.51 0.328 17.37 28.15 0.538 0.439 0.389 

0.1607 -15.92 1.000 235.5 30.00 33.2 -202.9 30.71 0.320 17.60 17.94 0.543 0.656 0.636 

0.1607 -15.88 1.000 235.0 28.00 6.7 -284.6 28.71 0.341 17.43 23.93 0.539 0.492 0.445 

0.1553 -15.99 0.911 370.3 29.35 35.3 -259.2 30.38 0.329 17.36 19.06 0.438 0.548 0.592 

0.1535 -15.99 0.794 313.3 30.35 57.2 -135.4 31.60 0.324 19.00 14.60 0.508 0.731 0.803 
Table Continued on Next Page. . . 
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Name 



Po 



ay m*/m K 



sym 

(fm-3) (MeV) (MeV) (MeV) (MeV) (MeV) 



(MeV) 



(MeV) 



(MeV) 



(fm) (fm) 



Ae-R (i° Ref. 
(fm) (fm) 



Skll 





1604 


-15 


95 





693 


242.8 


37 


52 


161.0 


234.7 


38 


20 





331 


17 


45 


11 


42 





541 


1 


126 


1 


253 





443 


[56] 


SkI2 





1575 


-15 


77 





685 


241.0 


33 


38 


104.3 


70.7 


34 


02 





316 


17 


03 


11 


68 





537 





951 


1 


096 





412 


[56] 


SkI3 





1577 


-15 


98 





577 


258.2 


34 


83 


100.5 


73.0 


35 


20 





293 


17 


77 


12 


77 





542 





908 


1 


045 





413 




SkI4 





1601 


-15 


94 





649 


248.0 


29 


50 


60.4 


-40.6 


30 


08 





344 


17 


48 


20 


83 





539 





665 





540 





582 


[56] 


SkI5 





1558 


-15 


85 





579 


255.8 


36 


64 


129.3 


159.5 


37 


01 





288 


17 


31 


10 


71 





536 


1 


034 


1 


316 





374 


[56] 


SkI6 





1591 


-15 


92 





640 


248.6 


30 


09 


59.7 


-47.3 


30 


64 





333 


17 


42 


19 


07 





539 





683 





602 





542 


[57] 


SLyO 





1603 


-15 


97 





698 


229.7 


31 


98 


47.1 


-116.3 


32 


67 





301 


18 


21 


16 


75 





580 





713 





727 





500 


[58] 


SLyl 





1603 


-15 


99 





698 


229.9 


31 


99 


47.1 


-116.6 


32 


68 





301 


18 


31 


16 


19 





582 





725 





752 





484 


[58] 


SLy2 





1605 


-15 


99 





698 


230.0 


32 


00 


47.5 


-115.2 


32 


69 





301 


18 


20 


17 


28 





579 





706 





705 





516 


[58] 


SLy230a 





1600 


-15 


99 





697 


229.9 


31 


99 


44.3 


-98.3 


32 


23 





301 


18 


20 


18 


19 





580 





668 





670 





525 


[59] 


SLy3 





1604 


-15 


97 





696 


229.9 


31 


99 


45.3 


-122.2 


32 


68 





300 


18 


21 


16 


74 





579 





706 





728 





490 


[58] 


SLy4 





1595 


-15 


97 





695 


230.0 


32 


00 


46.0 


-119.8 


32 


68 





301 


18 


24 


16 


60 





581 





712 





735 





490 


[60] 


SLy5 





1606 


-15 


98 





698 


230.0 


32 


01 


48.2 


-112.8 


32 


70 





302 


18 


28 


16 


06 





581 





728 





759 





481 


m 


SLy6 





1590 


-15 


92 





690 


229.9 


31 


96 


47.5 


-112.8 


32 


63 





300 


17 


53 


17 


12 





562 





689 





713 





487 


m 


SLy7 





1583 


-15 


90 





688 


229.8 


31 


99 


47.0 


-114.4 


32 


65 





299 


17 


35 


17 


45 





559 





680 





701 





490 


m 


SLy8 





1603 


-15 


97 





696 


229.9 


32 


00 


47.2 


-115.7 


32 


68 





301 


18 


18 


16 


40 





578 





717 





743 





485 


[58] 


SLy9 





1512 


-15 


79 





666 


229.9 


31 


98 


54.9 


-81.5 


32 


58 





300 


17 


59 


15 


93 





574 





742 





780 





486 


[58] 
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Name 



Po 



ay m*/m K 



sym 

(fm-3) (MeV) (MeV) (MeV) (MeV) (MeV) 



(MeV) 



(MeV) 



(MeV) 



(fm) (fm) 



Ae-R (i° Ref. 
(fm) (fm) 



SLylO 





1556 


-15 


90 





683 


229 


7 


31 


98 


38.8 


-142.3 


32 


62 





297 


17 


05 


20 


23 





549 





607 





608 





518 


[58j 


SkX 





1554 


-16 


05 





993 


271 


1 


31 


10 


33.2 


-252.2 


32 


36 





307 


16 


04 


19 


31 





476 





581 





620 





440 


m 


SkXce 





1553 


-15 


86 


1 


006 


268 


2 


30 


20 


33.7 


-238.5 


31 


46 





314 


16 


10 


18 


34 





480 





592 





634 





437 


[61] 


SkXm 





1589 


-16 


05 





966 


238 


2 


31 


20 


32.1 


-242.9 


32 


45 





306 


15 


69 


19 


56 





497 





597 





609 





469 




SkO 





1605 


-15 


83 





896 


223 


4 


31 


97 


79.1 


-43.2 


32 


86 





320 


17 


14 


14 


94 





546 





856 





815 





485 


[62j 


SkOp 





1602 


-15 


75 





896 


222 


4 


31 


95 


68.9 


-78.8 


32 


79 





309 


16 


07 


15 


74 





524 





802 





773 





453 


[62j 


SKRA 





1595 


-15 


78 





748 


217 





31 


32 


53.1 


-139.3 


32 


71 





308 


16 


71 


14 


79 





556 





766 





808 





455 


[63] 


MSkl 





1575 


-15 


83 


1 


000 


233 


8 


30 


00 


33.9 


-200.1 


30 


70 





318 


17 


21 


18 


08 





536 





651 





636 





492 


[64] 


MSk2 





1575 


-15 


83 


1 


050 


231 


7 


30 


00 


33.4 


-203.5 


30 


66 





318 


17 


22 


18 


06 





538 





654 





636 





494 


M 


MSk3 





1578 


-15 


82 


1 


001 


233 


6 


27 


99 


6.9 


-284.4 


28 


69 





339 


17 


08 


24 


37 





533 





485 





440 





562 


M 


MSk4 





1575 


-15 


79 


1 


050 


231 


2 


28 


00 


7.2 


-284.1 


28 


66 





339 


17 


06 


24 


04 





534 





494 





446 





563 


m 


MSk5 





1575 


-15 


79 


1 


050 


231 


2 


28 


00 


7.6 


-282.6 


28 


66 





339 


17 


07 


23 


90 





535 





496 





449 





561 


m 


MSkSs 





1561 


-15 


78 





800 


243 


8 


28 


00 


7.0 


-290.7 


29 


17 





337 


17 


25 


24 


15 





533 





481 





445 





549 


[65J 


MSk6 





1575 


-15 


79 


1 


050 


231 


2 


28 


00 


9.7 


-274.4 


28 


66 





339 


17 


09 


23 


10 





536 





510 





464 





556 


m 


MSk7 





1575 


-15 


79 


1 


050 


231 


3 


27 


95 


9.4 


-274.7 


28 


61 





340 


17 


11 


23 


12 





536 





509 





463 





557 


m 


MSk8 





1575 


-15 


79 


1 


100 


229 


4 


27 


93 


8.3 


-280.1 


28 


56 





340 


17 


12 


23 


38 





537 





507 





458 





561 


m 


MSk9 





1575 


-15 


80 


1 


000 


233 


4 


28 


00 


10.4 


-270.3 


28 


70 





339 


17 


15 


23 


02 





536 





510 





466 





553 


m 
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TABLE I - Continued 



Name 



Po 



ay m*/m K 



sym 

(fm-3) (MeV) (MeV) (MeV) (MeV) (MeV) 



(MeV) 



(MeV) 



(MeV) 



(fm) (fm) 



Ae-R (i° Ref. 
(fm) (fm) 



v070 

v075 

v080 

v090 

vlOO 

vl05 

vllO 

SKzl 

SKz2 

SKz3 

SKz4 

BSkl 

BSk2 

BSk2p 

BSk3 

BSk4 

BSk5 



0.1575 -15.80 1.050 231.3 28.00 -3.5 -361.7 29.49 0.337 17.41 24.61 0.541 0.466 0.436 0. 

0.1575 -15.80 1.050 231.3 28.00 -0.3 -342.0 29.32 0.337 17.31 24.40 0.539 0.477 0.439 0. 

0.1575 -15.79 1.050 231.2 28.00 2.3 -325.7 29.18 0.337 17.18 24.30 0.536 0.483 0.441 0, 

0.1575 -15.79 1.050 231.2 28.00 5.1 -304.3 28.94 0.338 17.08 24.22 0.535 0.489 0.443 0, 

0.1575 -15.79 1.050 231.2 28.00 8.8 -281.5 28.75 0.339 17.02 23.68 0.534 0.502 0.453 0, 

0.1575 -15.79 1.050 231.2 28.00 7.1 -284.6 28.66 0.339 17.05 24.11 0.534 0.492 0.445 0, 

0.1575 -15.79 1.050 231.2 28.00 7.5 -279.7 28.59 0.339 17.11 23.78 0.535 0.496 0.451 0. 

0.1600 -16.01 0.700 230.1 32.01 27.7 -242.5 33.60 0.294 17.35 19.03 0.554 0.616 0.641 0. 

0.1600 -16.00 0.700 230.1 32.01 16.8 -259.8 33.27 0.293 17.34 23.95 0.554 0.523 0.509 0. 

0.1600 -16.01 0.700 230.1 32.01 13.0 -242.0 32.77 0.293 17.35 28.05 0.554 0.465 0.435 0, 

0.1600 -16.01 0.700 230.1 32.01 5.8 -241.0 32.35 0.293 17.35 34.04 0.554 0.393 0.358 0, 

0.1573 -15.80 1.050 231.4 27.81 7.2 -281.9 28.47 0.342 17.22 23.72 0.539 0.499 0.449 0, 

0.1575 -15.79 1.042 233.7 28.00 8.0 -297.1 29.02 0.338 17.14 23.15 0.534 0.509 0.463 0, 

0.1575 -15.79 1.049 233.4 28.00 7.8 -298.1 29.02 0.338 17.11 23.20 0.534 0.508 0.462 0. 

0.1575 -15.80 1.123 234.8 27.94 6.8 -307.0 28.95 0.339 17.16 23.37 0.531 0.506 0.458 0. 

0.1575 -15.77 0.920 236.9 27.99 12.5 -266.0 28.90 0.339 16.91 23.01 0.530 0.511 0.466 0. 

0.1575 -15.80 0.920 237.2 28.70 21.4 -240.4 29.65 0.331 17.04 20.65 0.533 0.570 0.532 0, 
Table Continued on Next Page. . . 



547 [65 
552 [65 



555 [65 

561 [65 



560 [65 
563 [65 



559 [65 
490 [68 



549 [68 
596 [68 



643 [68 
566 [69 



552 [70] 

553 [701 



552 [71] 
551 [71] 
528 [71] 



TABLE I - Continued 



Name 



po CLV m*/m K 
(fm-3) (MeV) (MeV) 



(MeV) 



L 

(MeV) 



K 

(MeV) 



^(Po,l) 
(MeV) 



(MeV) 



(MeV) 



(fm) (fm) 



Ae-R (i° Ref. 
(fm) (fm) 



BSk6 





1575 


-15 


75 





800 


229 


2 


28 


00 


16.9 


-215.3 


28 


72 





339 


16 


74 


22 


73 





542 





511 





472 





564 


[71] 


BSk7 





1575 


-15 


76 





800 


229 


3 


28 


00 


18.0 


-209.4 


28 


69 





340 


16 


70 


22 


81 





541 





514 





470 





568 


[71J 


BSk8 





1589 


-15 


82 





800 


230 


3 


28 


00 


14.9 


-221.0 


28 


70 





341 


17 


12 


23 


13 





549 





506 





462 





578 


[72] 


BSk9 





1589 


-15 


92 





800 


231 


5 


30 


00 


39.9 


-145.4 


30 


62 





323 


17 


42 


17 


83 





556 





666 





643 





520 


[72] 


BSklO 





1593 


-15 


91 





920 


238 


9 


30 


00 


37.3 


-195.0 


31 


01 





321 


17 


51 


17 


48 





540 





671 





655 





491 


[73j 


BSkll 





1586 


-15 


86 





920 


238 


1 


30 


00 


38.4 


-189.9 


30 


99 





320 


17 


32 


17 


42 





536 





671 





658 





486 


[73j 


BSkl2 





1586 


-15 


86 





920 


238 


1 


30 


00 


38.0 


-191.4 


31 


00 





320 


17 


30 


17 


51 





535 





668 





655 





487 


[73] 


BSkl3 





1586 


-15 


86 





920 


238 


1 


30 


00 


38.8 


-188.0 


30 


99 





320 


17 


34 


17 


34 





536 





673 





661 





485 


[73] 


BSkl4 





1586 


-15 


85 





800 


239 


4 


30 


00 


43.9 


-152.1 


30 


92 





323 


17 


17 


17 


22 





537 





680 





666 





492 


M 


SK255 





1573 


-16 


33 





797 


255 





37 


40 


95.1 


-58.3 


38 


77 





268 


19 


28 


12 


44 





561 





967 


1 


153 





389 


[75] 


SK272 





1553 


-16 


28 





773 


271 


5 


37 


40 


91.7 


-67.8 


38 


71 





264 


19 


21 


13 


48 





545 





916 


1 


068 





395 


[75] 


QMCl 





1373 


-14 


00 





926 


328 


7 


29 


68 


-6.7 


-504.2 


31 


03 





272 


23 


72 


21 


44 





663 





708 





555 





680 


[76] 


QMC2 





1403 


-14 


29 





834 


330 


1 


28 


70 


8.7 


-408.4 


29 


84 





290 


18 


84 


21 


00 





526 





586 





544 





520 


[76] 


QMC3 





1607 


-15 


98 





825 


366 


9 


45 


78 


91.8 


-211.0 


47 


07 





197 


21 


17 


18 


92 





508 





847 





921 





418 


[76] 


KDEOv 





1608 


-16 


10 





717 


228 


8 


32 


98 


45.2 


-144.9 


33 


80 





291 


17 


27 


18 


22 





552 





672 





689 





488 


[72 


KDEOvl 





1646 


-16 


23 





744 


227 


6 


34 


58 


54.7 


-127.2 


35 


34 





281 


17 


66 


17 


25 





558 





730 





757 





474 


im 


LNS 





1746 


-15 


31 





826 


210 


8 


33 


43 


61.5 


-127.4 


34 


63 





276 


15 


77 


14 


10 





517 





765 





878 





404 


m 



D. Local Approximation 



We now consider the Hohenberg-Kohn (HK) functional in the context of uniform and 
near-uniform matter. For a uniform p, given the requirement of translation invariance, we 
may notice that S in Eq. (1201) must be a symmetric function of the difference of its arguments, 

S{p,r,ri) = S{p,r-ri) = S{p,ri - r) . (70) 

Upon comparing (1201) and (1551) . we arrive at the constraint in uniform matter 

drS{p,r)= JdrS{p,r-ri) = ^. (71) 

Given the evidence for a local nature of nuclear energy functional in the energy formula and 
the above constraint, it should be possible to approximate the operator S for uniform p by 

5(p,r)^^5(r), (72) 

in the action of this operator on sufficiently uniform densities pnp- Also, such an approx- 
imation should be valid for densities p that are not uniform generally, but are sufficiently 
uniform locally, with the operator then given by 

S{p,r,r,)^^^6{r-r,). (73) 

On applying an inverse operator to S for uniform p, to the middle and r.h.s. expressions 
in (I7T|) . we find a constraint on the inverse operator, 

I drS-\p,r) = ^^. (74) 

A consistent approximation for the inverse operator, in matter sufficiently uniform locally, 
is 

S-\p,r,r^)c^^^^5ir-r,). (75) 

From ( H9l) and ( 17^ . it follows that the isovector density in uniform matter can be repre- 
sented as 

Pa = ^^P + Oiv'). (76) 

The isovector density is then seen normalized locally in such a manner that it approaches 
the normal density po when isoscalar density itself approaches po. Equation ( l76l) is further 
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expected to hold in nonuniform matter that is sufficiently uniform on the scale of nonlocality 
in the energy functional. The contributions from such regions to A/aa{A) in flSTl) are then 

^ ^dr-i^ + ..., (77) 



aM) J Sipir)) 

where dots represent contributions from other regions. In (1771) . the contributions from locally 
uniform regions combine as capacitances of independent capacitors in electrostatics, with 
p/2S{p) playing a role of local density of capacitance and A/2aa{A) playing the role of the 
net capacitance, cf. Eqs. ( l29l) and ( l30l) . 

One extreme limit that helps to develop an intuition, following the local approximation, 
is that of a constant symmetry energy, S{p) = a^. We find in that limit, as common sense 
suggests, that aa{A) = , cf. f l77|) . In the context of Eq. f|T4l) . we can see that aa{A) and 
are expected to coincide when af approaches infinity. As a further implication of 5* = a^, 
we find, cf. (1391) . that the net nuclear density becomes identical to the density for symmetric 
matter up to the 4**^ order in asymmetry not just 2"^^, p(r) = po{r)+0{r]'^). For the isovector 
density, cf. (1761) . we find p°(r) = poir) and Pair) = p{r) + 0{ri'^) = Po{r) + 0{rf). Since 
the densities coincide for symmetric matter, p°(r) = poir), the condition fl53l) . signaling 
a breakdown of the quadratic approximation to the symmetry energy, is never satisfied. 

These S = results suggest that, for interactions for which the function S varies 
weakly with p, the isoscalar density should exhibit a weaker variation with asymmetry, 
than the isovector density. The result ( 1761) suggests that the stronger variation of S{p) 
with p, the more pronounced should be the differences between the isovector and isoscalar 
densities. Moreover, Eq. (177j) suggests that the stronger that variation the larger should be 
the difference between aa{A) and . If that difference were attributed to the nuclear surface, 
a larger difference would correspond to a reduced value of surface symmetry coefficient af , 
cf. (1771) . Finally, the stronger the variation of S with density, the sooner the invariance of 
the isoscalar density, and then of isovector density as well, should break down. Indeed, in 
a region where the local approximation holds, we find, from (1391) . 



, 2 r A\ 

Po + r] aa{A) 



Sp{r 



A daa ^ agjA) / _ po{r) dS 



0{v') ■ (78) 



.a„(A)dA S{po{r))\ S(po(r))dp, 

From (ITS]) , it appears that, simultaneously, a deviation of the local symmetry coefficient S 
from the global coefficient Oa, local dependence of S" on p and the dependence of on A 
(itself induced by the dependence of 5* on p) , all conspire to produce a significant coefficient 
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of quadratic expansion for the isoscalar density around a given location r. Further evidence, 
for a reduced range of vahdity of the imphcations of the second-order ii^a-expansion, for a 
stronger p-variation of S, is the fact that the enhanced differences between the isovector and 
isoscalar densities should lead to an earlier breakdown of the condition flS51) . whether when 
staying at a given location in the matter and increasing rjy or when departing from nuclear 
surface at a fixed rjv- 

While the po-parameter L, Eq. fl62|) . is commonly used to quantify the p-dependence 
of S, obviously variation of S over a wider range of subnormal densities must impact the 
properties of matter associated with the asymmetry. In the local expressions ( 1761) and ( 1771) . 
it is seen that the inverse S is weighted with density. In the practical calculations of half- 
infinite nuclear matter of the next section, we find that local values of 5* impact the isovector 
density and aa{A) within the isoscalar density range of po/4 < p < po- For positions in the 
self-bound matter where p < Po/4, the operator S becomes delocalized and, thus, any 
potential impact of the specific values of S{p) gets smeared out. Since any realistic S{p) 
eventually tends to zero as p ^ 0, a negative value of L at po, see Fig. Ej leads in practice 
to an S with a weaker variation with density over the pertinent region, than a similar in 
magnitude but positive L. 

Selected results, of the local approximation pursued here, have been arrived at in the 



past following the local 
Usmani 



16 



24 



79| . Specifically, Bodmer and 



homas-Fermi approximation 
Tyl (see also 2M) have demonstrated, in the absence of Coulomb effects, an inverse 

n n 

proportionahty of pnp to the symmetry energy at a local net density. In [24] (see also [16[ 
by Steiner et ai), the ratio A/aa{A), i.e. up to a factor the capacitance of a nuclear system 
for asymmetry, was found to be given by the integral on the r.h.s. of Eq. (177|) . 



E. Functional Continuity and Kohn-Sham Methodology 

In employing the Hohenberg-Kohn functional, we found it important to invoke a smooth- 
ing and interpolation procedure making the resulting functional analytic in nucleonic densi- 
ties. If one were to find the Hohenberg-Kohn functional in reality, though, e.g. progressing 
over certain directions in the density space, the functional would be defined for discrete par- 
ticle numbers only and would, otherwise, exhibit definite discontinuities across the density- 
ranges of interest. The latter are implied by the abrupt changes in experimental ground-state 
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energies, superimposed on top of a smooth variation of the energies with nucleon numbers. 
Those abrupt changes can be, in particular, attributed to filhng of single-particle energy 
shells. In Kohn-Sham methodology 80|, the difference between the actual (in most prac- 
tice hypothetical) Hohenberg-Kohn functional, on one hand, and the kinetic energy of such 
form as for a noninteracting system, on the other, is represented by a continuous functional 
of density. The advantage of the Kohn-Sham methodology is that the assumed functional 
yields the exact theory in the limit of no interactions. Otherwise, discontinuities in the ap- 
proximate functional, due to its kinetic part, on top of what would follow from using discrete 
particle numbers, may bear similarities to those for the exact Hohenberg-Kohn functional. 

The Skyrme-Hartree-Fock model (SHFM), discussed in the next subsection and well-tried 
in nuclear physics, see e.g. j2| and jsil], expands the Kohn-Sham methodology, in that the 
kinetic-energy terms for nucleons in the approximated functional get multiplied by factors 
dependent on density. The latter allows for some control over the discontinuities resulting 
from the functional. The remaining part of the functional is kept only minimally nonlocal 
through the use of lowest-order terms in density gradients. SHFM stems from modeling 
internucleon forces; a strategy relying from the outset on a Kohn-Sham type functional has 



82|. 



been pioneered in nuclear physics by Fayans et al. 

In an application of the SHFM, later in the paper, to the domain of semi-infinite matter 
at changing asymmetry, it might seem that no smoothing of the approximate energy func- 
tional is required to make our preceding considerations, based on a functional expandable 
in asymmetry (cf. Eq. ( 120|) ). directly applicable. Such a conclusion might be inferred from 
the fact that the half-infinite matter is characterized by infinite nucleon numbers and by 
a continuum energy spectrum. However, it will turn out that, on account of the so-called 
Friedel oscillations extending from nuclear surface into the matter, the r^-derivatives of nu- 
cleon density from functional minimization, of order third and higher, steadily increase with 
distance away from the surface. The implication of this behavior is that the expansion in 
asymmetry, which we have considered for the functional, must be just asymptotic in the 
specific domain. A direct application of our Hohenberg-Kohn considerations to SHFM for 
half-infinite nuclear matter, without any smoothing, could be then regarded particularly 
interesting as an extrapolation of our considerations. On the other hand, regarding the 
formal issue of smoothing and of functional analyticity, any smoothing that would lead to 
a smoothing of the resulting densities over a small range of Fermi wavevectors, would tame 
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the mentioned derivatives with respect to asymmetry. In contrast to half-infinite matter, 
the SHFM functional apphed to uniform matter renders fully analytic results for energy 
without any smoothing. 



F. Skyrme-Hartree-Fock Model 



The Skyrme-Hartree-Fock model gets exploited in this paper in different ways. On one 
hand, the model has been already utilized to illustrate properties of uniform matter. On 
the other hand, the model will be utilized in the half-infinite matter calculations. Results 
of those calculations will be exploited to verify assertions concerning the Hohenberg-Kohn 
energy functional. Since the SHFM energy functional is expressed in terms of single-particle 
wavefunctions rather than in terms of densities and Pp, the inferences regarding 

the HK functional can be nontrivial. 

In SHFM, the energy is 

E — -Bskyrme + Ec + -Bpair , (79) 

where the two last terms Eq and i^'pair are, respectively, the Coulomb and pairing energies. 
The Skyrme functional is 



-Eskyrme = / dr Cskyrme (j") 



+ 



12 



+ 



+ 



t2 

4 

3^1 
16 



i + f) (vp)'-(i+-.)E(vp,)" 



1 , 



(80) 



-hpVJ-b'^ ^P.VJq 
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where the g-index summations are over neutrons and protons, t = Tn + Tp and J = Jn + Jp- 
Particle densities and scaled kinetic energy and spin densities are obtained from summations 
over the single-particle wavef unctions: 



a 

r,ir) = J2KV<Pi{r)VMr) 

a 



(81) 
(82) 
(83) 



Here, {n^,} are occupation numbers of the single-particle states. In less elaborate Skyrme 
parameterizations, the two spin-parameters are combined into one: 



h = K = U/2 . 



(84) 



The gradient terms in the interaction part of the energy functional (l80l) (multiplied by 
adjustable constants), represent a simplified form of interaction nonlocalities. 

Minimization of the SHFM functional (179|) . with respect to the wavefunctions, produces 
self-consistent equations for the latter 



iK + Uc){r) 



(85) 



In the equations, Uc is Coulomb potential and the nuclear single-particle hamiltonian hg is 
of the form 



Kir) 



-V Ba(r) V + Uair) - iWair) (V x a) 



(86) 



The mass potential in hg is given by 



5„ 



2m* 2m 4 



+ 



4 



(87) 



and m* is position-dependent effective mass. The potential Ug in fl86|) is 



+ 



12 



(2 + a)(l + f)/ 



+ ^3) (a XI ^9' +2pPg 



1 , 



(88) 



l^+X2]Tg 
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Finally, the form factor Wq in (l86i) is 

Wg = hVp + b',Vpq. (89) 

Above, we suppress the terms in the SHFM functional and single-particle hamiltonian that 
are normally small, such as those proportional to the square of spin density. 

In the context of investigating the effects of neutron-proton asymmetry, we shall introduce 
average and deviation-from-average quantities, for the terms in the single-particle hamilto- 
nian and for the hamiltonian itself. These will be denoted, respectively, with either the lack 
of species index or with that index replaced by 'a'. (Unlike in the case of density, no special 
normalization will be adopted here for the quantities with the index 'a'.) Thus, e.g. for the 
mass potential, we will have 

B = ^J1±3l = ^ + A \3t, + (5 + 4x2)^2] , (90) 
2 2m 16 L ^ ' ^ ^ 



and 



Ba = ^^V^ = IT + 2^2)t2 - (1 + 2xi)ti] . (91) 
/ 16 



For the effective mass, we will have 



where the approximation holds for low asymmetries, and 

1/1 1 \ B„ m*B. 



"'^"4Un 2{B^-Bl)- B ■ ^^^^ 



For the optical potential, we will have, one hand, U = {Un + Up)/ 2 and, on the other, 

(94) 



Ua =^1^—Ep = _^ (1 + 2Xo)pnp - ^ (1 + 2a;3)p'' Pnp 



1 h' 

+ — [(1 + 2X2)t2 - (1 + 2Xi)ti] Tnp - -j- VJnp , 



where r„p = r„ - Xp and J„p = J„ - Jp. 

In uniform matter, the kinetic energy densities are 

3 



5 ' 

with the Fermi wavenumbers given by 

,1/3 



TPqk', (95) 



kFq={S7T^P<j)' . (96) 
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In symmetric matter, the partial densities are equal and, thus. 



Pq = p/2 and Tg = r/2 . 



(97) 



5 /3 

With Tg oc pg , the partial kinetic energies may be expanded in asymmetric matter around 



their symmetric value. 



1 ± 10 ( pnp 

3 2p 9 V 2p 



± . . . 



(98) 



where kp is wavevector for the symmetric matter at density p. Upon inserting the symmetric 
values into (1801) . the energy per nucleon of symmetric uniform matter becomes 

(99) 



En , , "if? , o 3to , ^3 ^-Li 3 



For a future use, from (pHll . the difference in partial kinetic energies may be represented at 
small asymmetries as 

Tnp^Pnpkl. (100) 

Finally, upon exploiting (p8|) in the energy (IHOl) for uniform matter, the density-dependent 
expansion coefficient in asymmetry becomes 



6m 



kl-^j{l + 2xo)p-^{l + 2xs)p 



a+l 



1 

24 



[3tixi - {A + 5x2)t2]pkl. 



(101) 



As is apparent, for uniform matter the same powers of density appear in the symmetry 
energy S and the energy per nucleon Eq/A of symmetric matter. From the interaction con- 
stants, the dimensionless constants xq, xi and 0:3 affect the strength of the three interaction- 
terms in S, with different density powers, without affecting Eq/A. From among the different 
SHFM quantities, the mass potential has a relatively simple linear density dependence. As 
should be obvious from Table [H the number of Skyrme interaction constants in use in the lit- 
erature is large and it is, obviously, growing. For the values of constants, we refer the reader 
to the references indicated in the table. Otherwise, one larger compilation of constant values 
can be found in Ref. 32 1. 



A combination of SHFM with the semiclassical expansion for kinetic energy 
of density. 



in terms 



(102) 
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valid for slowly varying densities, allows one to gain an insight into the symmetry opera- 
tor S and its inverse. For simplicity, we assume the spin densities Jq to be zero. Upon 
inserting (11021) into (Mil) and employing the identity 



(103) 

Pnp=0 



we find for the operator 

S{p,ri,r2) ^So{p,ri)6{ri -rs) 

dr Si (p, r) — • [6{r - r i) V5(r - rs) + 6{r - ra) V5(r - ri)) (^194) 

+ j drS2{p,r)VS{r-ri)-V5{r-r2), 



where 



So = - + -^(^X {- + l[{2 + x,)t, + {2 + X2)t2]p] 
p 18p \ p / [ m 4 '- -'J 

'5i = -^{^ + ^ [(2 + xi) ti + (2 + X2) t2] p} , 



(105) 



Those nonlocal gradient terms in the combination of results (11041) and (llOSp . which are 
multiplied by force constants, are associated with the nonlocality of interactions. On the 
other hand, those gradient terms, which are multiplied by h"^, are associated with quantum 
nonlocality. It may be noted that, for an approximately uniform p, the operator S given 
by (I104p becomes a function of (ri — ^2) only. In addition, when both p and p„p onto which S 
acts approach uniformity, then S given by (I104p approaches the local result (!73l) . 



IV. HARTREE-FOCK CALCULATIONS FOR HALF-INFINITE MATTER 
A. Skyrme-Hartree-Fock Model for Half-Infinite Matter 

With the interest in low-density features of symmetry energy, we now turn to an examina- 
tion of half-infinite nuclear matter {4, gI within SHFM. The half-inifinite matter is uniform in 
two cartesian directions and generally nonuniform in the third that we choose to be z. Given 
the transverse uniformity, the solutions to the SHFM equations (!85|) may be looked for as the 
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eigenstates of transverse momentum, 0(r) = "^{z) exp{ik±r±). Here, k± is the transverse 
wavevector and \I/ is a spinor depending on z. In the single-particle hamiltonian fl86|) . only 
the last-r.h.s. spin-orbit term has a vector character. With the form factor Wg of Eq. (15^ 
pointing in the 5;-direction, that spin-orbit term, acting on a wavefunction, becomes 

- iWy{r) (V X a) 0(r) = -tWg{z) (V x a)^ 0(r) = Wq{z) {k± x a), 0(r) , (106) 

where 

iy,(z)=64^ + &l^. (107) 
dz dz 

If we next orient the x-axis along the A;_L-direction, see Fig. HJ we find that the spin-orbit 
term in the hamiltonian becomes proportional to the Pauli matrix along the y-direction: 

With this, it turns out that the eigenstates of the hamiltonian may be searched for as the 
eigenstates of operator of spin projection onto the y-direction, recognized as the direction 
of orbital angular momentum calculated relative to the far-away interior of matter: 

Mr)=i^k^kM^)Xxe'''^''^. (109) 

Here, xx is an eigenstate of ay, 

o-s,Xa = Axa, (110) 

the eigenvalues are A = ±1, and ip{z) is a scalar wavefunction. 

From ( l85l) and f ll08l) -( fTT0|) . the wavefuctions satisfy the differential equation 

e.,^...TMi^)=[-^^^^^^ + ^^ (111) 

The single-particle energy above can be represented as 

'k, = ir-^ + U,oo, (112) 

^"'(JOO 

where the quantities with the infinity index refer to a uniform interior of the matter far away 
from the surface, and where the wavevector squared is 

e = ki + {k^y . (113) 

The last two equations principally define k'^ > in terms of e. The 2;-component of the 
wavevector must be defined in the asymptotic region in the matter, as the z-component of 

39 



0.4 r 




0.35 - 



p (fm-3) 



0.15 ^ 



0.05 - 



0.25 - 



0.1 - 



0.2 - 



0.3 - 












z(fm) 



8 



9 



FIG. 4: Within our half-infinite matter SHFM-calculations, the z-axis is ahgned with the direction 
of nonuniformity of matter, while the x-axis is aligned with the single-particle transverse momen- 
tmii. The single-particle particle states are chosen to have a definite projection of spin onto the 
remaining y-axis. 

single-particle momentum is generally not conserved, due to matter nonuniformity. Other- 
wise, in the asymptotic region, the eigenstates of single-particle hamiltonian combine waves 
moving towards and away from the boundary of matter, following a reflection from that 
boundary. 



Here, the individual plane waves in the combination are normalized to unity and 6 is the 
phase of the incident relative to the reflected wave. The phase shift depends on the choice of 
origin for the z-axis. However, changes in the phase shift with species or with momentum, 
such as in 6, do not depend on the choice of origin. The arbitrary overall phase-factor 
multiplying the wavefunction in flll4p is chosen so as to make the wavefunction ip real 
asymptotically. With the wavefunction and its derivative being real, and with the coefficients 
in the wave-equation (11 111) being real, the wavefunction ip, after setting the phase-factor, 
becomes real everywhere. Far outside of the nuclear surface, in the classically forbidden 
region, the wavefunction must decrease exponentially with distance away from the surface. 




(114) 
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where we assume that the matter is positioned towards the negative direction of the 2;-axis, 
and 

^k,= ^kl-^-^. (116) 

Provided self-consistent solutions for the half-infinite matter may be found, different 
observables can be obtained by carrying out integrations over wavefunction continuum, with 
the summations over discrete states, such as in (jSTl) . specifically replaced by 

-1 /» /"OO 1 /*1 /"OO 

E(-) - E d*r (.) = E /„ dcos« I (.) . (117) 

ct A A 

in accordance with the normalization of the wavefunctions in flll4p . On the r.h.s., the 
wavevector components cos 6 and kj_ = k sin 9. In particular, with the Hartree- 

Fock occupations of nl = 9{kp — k), the nucleon and kinetic-energy densities take, respec- 
tively, the forms 

Pii^) = y27T^ dcos^ / dkk"^ \tlJk^k^Xq{z)\'^ , (118) 



and 



^g(^) = X^^^^^ dcosO J dkk^ \ k^ | V'fcxfcf Ag(^) T + 



k'' / 1 ; 2^ 

,2 ( l2 Li. /^\^ 



dz 



(119) 



Given the symmetry of half-infinite matter, only the 2;-component of spin density can be 
finite. With (IH3]), ffT09D and ffTTOll . cf. Fig. H that component takes the form 

^'('^) = 5^7^ / [ ^dkk^k^ \i^k^k^xq{z)\\ (120) 



B. Numerical Solution of the SHFM Equations 

As the first step, we solve the SHFM problem for nuclear matter for a certain assumed 
asymmetry rjy in the interior of the matter. Specifically, under the constraint of fixed riv, 
with = py {1 ±riy)/2 and employing fl95|) and fl96l) . we seek pv that minimizes E/A of 
Eq. (IHOj) for uniform matter. That solution provides features of the semi-infinite matter in 
the asymptotic region. 

Given the densities p^^ and the corresponding Fermi momenta k^^, we next discretize the 
wavevector space. Specifically, aiming at a calculation of the integrals f lll8p -( |T20l) for observ- 
ables, we discretize the space according to the abscissas {uf}fLi for Gaussian integration 
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in angle, 

cos9j = ^{l + Uj), j = l,---,M, (121) 
and according to the abscissas {tf^}^^^ for Gaussian integration in wavenumber, 

ihgf = \{l + U) {kp.f , z = 1, ■ ■ ■ , iV . (122) 

The choices fll2ip and (11221) for the discretization are made to provide a uniform coverage 
of wavevector space, when integrations such as in (I118p - (I120I) are approximated numerically 
by sums, following 

jocose l\kk^.) = ^f jjuAt{.)^^^Y.^fwf{.). (123) 

Here, w are the respective Gaussian weights. The wavevector components for the expressions 
summed over on the r.h.s. of (11231) are then 

ie^g = h,cosdj, (124) 

and 

kl^ = hqsine,. (125) 

The wavefunction equations (II 111) are solved self-consistently for the chosen discrete wavevec- 
tor values, at two spin orientations. 

We integrate the differential equations (lllip . for ip at individual wavevectors, by start- 
ing out far out from the surface. Over there, we relate the derivative to wavefunction in 
accordance with the asymptotic relation (IllSp . i.e. impose the condition 

^\z) = -Ki){z) . (126) 

and we move with our solution towards the matter and into its interior. Formal solutions 
to the differential equations (1 11 II) generally consist of two components: 

ij{z)=i)-{z)+i)^{z). (127) 

In the above, the ip' component drops in magnitude when moving away from the nuclear 
surface, according to the asymptotics of Eq. (11151) . That component represents, up to 
a factor, the solution we are looking for. On the other hand, the ip^ component grows when 
moving away from the surface, behaving as 

V'+(^) oc e+'^^ (128) 
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Have we moved away from the surface in integration, any inaccuracy in fll26p compared to an 
exact relation for ip' , or any accuracies in the integration, would admix the component 
to our numerical solution. The component would eventually begin to dominate over ip~ 
due to its exponential growth. However, if we move towards the matter, the opposite takes 
place, i.e. any admixture fades away due to its exponential decrease. For such reasons, 
in fact, it is not very important to insist on fll26p in the starting conditions, as long as any 
results very close to the starting point are not utilized. E.g. one could start integration with 
ip' = 0. Due to the linearity of the wavefunction equation, further, the starting magnitude 
of the wavefunction is not even important, i.e. e.g. ip = 1 could be principally used. This 
is because the wavefunction normalization can be reset after the integration arrives into the 
interior of matter, by imposing the condition 

2 



1 d^k 



(129) 



kf dz 

following (11141) . In practice, because of the Friedel oscillations that will be discussed, which 
die off slowly in z, we prefer to impose the condition (11291) with the l.h.s. in the condition 
averaged over a distance in z of ~ rc/kp. 

An attractive high-accuracy method for solving second-order differential equations, such 
as the Schrodinger equation, of fifth order in the employed step, is the Cowell-Numerov (CN) 



method 



84,^ 



86l |. High order of a method generally allows for large steps in integration 
preventing accumulation of rounding errors, when numerous small steps need to be taken. 
The CN method applies directly, though, to the equations where the second derivative is 
expressed exclusively in terms of the independent variable and the integrated function. For 
those Skyrme interactions in Table [U for which m*/m ^ 1, the position-dependent effective 
mass in the wavefunction equation (II 111) introduces, however, the first-order derivative of 
wavefunction into the single-particle equation, preventing a direct application of the CN 
method. With Bq(z) = h^2mq{z), we specifically get for the derivative term in (lllll) : 

±B-^P = B^ + — ^ (130) 
dz dz dz^ dz dz 



A straightforward work-around has been put forward by Dobaczewski et al. 871]. Thus, on 
introducing 

Vk^kf\q{z) = By^{z) 'ipk^k^Xgiz) , (131) 
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the r.h.s. of (11301) can be rewritten as 



51/2 ^ + 



where 



Upon rewriting of the derivative term in (11321) . the first-order derivative disappears from the 
expression at the cost of the appearance of a term hnear in the function, but the latter is not 
a problem for the CN method. Upon inserting the derivative term into (lllll) . the differential 
equation for v takes the form 



where 



Uf{z)=U,{z) + /\U,{z). (135) 



The equation for v above can now be solved with the CN method and ip can be obtain 
from V following (11311) . The asymptotic z — » ±00 behaviors are of the same type for v as 
for if), as B approaches a constant either inside or outside of the matter. 

In the discretization of wavevector space, we typically utilize abscissa numbers of M ~ 14 
and ~ 28. The typical spatial step that we employ, when integrating the wavefunction 
equations with the CN method, is /S.z = 0.02 fm. We start from a distance Zmax displaced 
from the matter surface position of Zq by z^^^^ — ;zo — 14 fm and we continue the integration 
of the wavefunction equations down to a position Zmin in the matter, away from the surface 
by Zq — z^i^ ~ 16 fm. After each solution pass, the local particle and kinetic energy densities, 
of Eqs. (IllSp and (11191) . are updated, with the integrals (I117P in the densities represented 
in terms of the summations ( 1123^ . 

For initialization of the solution iterations, starting densities must be assumed and we 
employ the combination of Fermi shapes: 

Pn,pW = J[p^W+<;(^)], (136) 



where 



1 + exp|_(2; - zoj/d\ 
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and 

Pnpi^) = sgn(prp)min(|p^p(2;)|,p(z)) , (138) 

with 

= r/"~''\^ W.1 • (139) 

1 + exp [{z - Zq- ARa)/da\ 
In the above, zq is the starting position of nuclear surface, d is the diffuseness of isoscalar 
density, ARa is the displacement of isovector density relative to the isoscalar density, and da 
is the diffuseness of isovector density. The role of Eq. fll38p is to ensure that the magnitude 
of asymmetry does not exceed 1 locally. The initial kinetic energy densities for the wave 
equations are estimated from the densities using Eq. (11021) . The spin densities are at first 
assumed to vanish 6| and are thereafter found to rise to finite values during iterations. 
The surface diffusness is first taken in the initialization equal to d ~ 0.55 fm. After self- 
consistency is first reached for symmetric matter, the value of d is read off from the obtained 
p{z) and that value is then used to initialize a restarted set of iterations for the given 
interaction. For calculations of asymmetric matter, the initial values used are ARa = 0.6 fm 
(cf. Table [T]) and da = d. After the self-consistency is first reached for a small asymmetry r]v, 
the values of ARa and da are read off from the small-T^y solution and used subsequently in 
the initialization of all other calculations of asymmetric matter for the specific interaction. 

Within iterations for self-consistency, an attempt to use directly output densities from 
an J'th solution iteration, [p]^^^ and [t]^^^, as input densities to the (J+l)'th iteration, 
generally introduces problems of two types. One problem is an undesired gradual shift of the 
solution in the direction outside of matter, towards a boundary of the computational region. 
Another problem are instabilities. To remedy both problems while optimizing convergence, 
as an input to the subsequent iteration, we employ, see also , a linear combination of the 
preceding input and of the shifted output, 

[p,(z)][;^ = f3{z) [p,{z + u,)]l^ + [1 - f3{z)] [p,{z)]l , (140) 

with the same P{z) and uj for pg and r^. In the practice of setting the value of uj, we 
found the requirements of either an unchanged nucleon number within the computational 
region, or of a fixed position for the net density reaching half of its asymptotic value, to 
be comparably effective in speeding up the iteration convergence and in preventing the 
undesired drift. The factor (3 has been made to vary smoothly across the computational 
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region, and to reach larger values, up to 0.4, outside of the matter and smaller inside the 



matter. For two of the Skyrme parameterizations, SkzO and Skzml [68|] , we could not arrive 
at consistency at finite asymmetry, no matter what f3{z) we employed. For any of the 
Skyrme parametrizations, the use of an excessive zq — Zj^m might induce long- wavelength 
instabilities in the solution, inside the matter, evidenced by a convergence of the solution 
iterations, but to slightly different density shapes Pn,p{z) in the matter, depending on the step 
A2; and on the starting position Zmin- Except for the immediate vicinity of the boundaries 
of the computational region, the solutions for interactions represented in Table [H for our 
choices of -Zmin.max and A^, are stable with respect to variations in the boundaries and in 
step magnitude. When comparing to each other our solutions for different interactions and 
different asymmetries in the next section, we reposition the solutions, so that the location 
where the net density drops to half of its asymptotic value does not change from one solution 
to another, p{zq) = /9y/2. 



C. Role of the Spin-Orbit Coupling 

Figure [5] shows the typical central potential U, spin-orbit form factor Wz and net density p, 
arrived at in the course of solving the SHFM equations until consistency, in symmetric half- 
infinite nuclear matter. As seen in the figure, the density p and potential U exhibit different 
diffusenesses, with p being steeper. Otherwise, consistently with ( |89l) and with general 
physical requirements, the spin-orbit form factor is finite only in the surface region. 
With this, it can affect macroscopic features of the nuclear surface. 

Quantitatively, with k± kp ^ 1.33 fm~^, the spin-orbit term k_LWz is small in the 
wavefunction equation (11 111) , lower by as much as a factor of ~ 4 or more at the maximum 
of its magnitude compared to the central potential U, see Fig. El The sign of Wz makes the 
spin-orbit potential attractive for spins parallel to the orbital angular momentum, A = -|-1, 
and repulsive for antiparallel, A = — 1. On account of the potential, nucleons with A = +1 
get, on the average, displaced outward relative to the nucleons with A = — 1, as illustrated 
in terms of the respective densities in Fig. [51 This displacement, in turn, produces a positive 
spin density Jz in the surface region, cf. ( I120p . A finite spin density modifies the central 
potential which involves a divergence of that density, see Eq. fISS]) . With 64 + > for 
a Skyrme parameterization, the correction to the central potential is attractive in the more 



46 




z (fm) 

FIG. 5: Results of reaching self-consistency in symmetric half- infinite matter for the SHFM equa- 
tions with the SkMP interaction, shown as a function of position along the axis perpendicular to 
the surface. The top panel displays the central potential U and the spin-orbit form factor Wz- 
The bottom panel displays the net nucleon density po and partial densities of nucleons with spin 
parallel Pq^^ and antiparallel Pq~~^ to orbital momentum. Also displayed is the net density Pg^^"'^ 
obtained when the spin-orbit potential is put to zero, Wz = 0. 

inner part of the surface and repulsive in the outer. This steepens the central potential and, 
in consequence, makes the net density steeper in the surface region compared to the case 
when the spin-orbit term is ignored, see Fig. [5] where the density arrived at without the 
spin-orbit term is also shown. 

A steeper profile of net density due to the spin-orbit term has been first predicted by 
Stocker, within a variational consideration [88;]. Details behind his prediction, however, differ 
from our results. Thus, he has anticipated a steeper density-profile for A = —1 than A = 1 
nucleons, while we find the opposite to hold, see Fig. [51 We next turn to the systematics of 
density with changing asymmetry. We will come back to the role of spin-orbit coupling in 
the surface region in the context of coefficients characterizing surface energy. 
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D. Densities for Different Asymmetries 

In Section UTTl we have indicated that, in consequence of the quadratic dependence of Ea 
on pnp{r), the normahzed difference of nucleon densities for a given mass number, pa{A){f), 
as well as the net density, p{r), should both weakly depend on asymmetry. Figures [H] and [7] 
show, respectively, the isoscalar and isovector densities from solving the SHFM equations for 
half-infinite matter at several rj, for sample Skyrme parameterizations. It is apparent in the 
figures that the anticipated scaling is generally very well satisfied. The scaling behavior holds 
even at r/ = 0.3 for the MSk9 and SkMP interactions, when the drip asymmetries for those 
interactions are close-by, at riy = 0.339 and 0.336, respectively. For the SkI5 interaction, 
the asymmetry rj = 0.3 is already above the drip asymmetry of 0.288 and the corresponding 
densities are not displayed in the figures. A closer examination of the densities in Figs. [6] 
and [7] reveals that, in the interior of matter, the scaling worsens a bit with an increase in the 
magnitude of L and, in the density tails, the scaling worsens with an increase in the value 
of L. These findings are consistent with the expectations developed in Subsection IIIIDI 
Regarding implications of Eq. fl75]) . the overall symmetry coefficient for half- infinite matter 
is identical to that for the volume, = a^. Consistently either then with (1751) or with (IMIl . 
in the interior of the matter in Fig. [6l the net density goes up or down depending on the sign 
of L, to the order of ■i]'^. Further, consistently with the expectations developed in the context 
of the local approximation in Subsection IIII Dt in the density tails, the scaling is seen to be 
a bit worse satisfied for the isovector densities in Fig. [TJ than for the isoscalar densities in 
Fig. O Nonetheless, one can claim that, in the end, the scaling for the isovector densities is 
more impressive than for isoscalar. This is because the isoscalar densities are rather bland 
across the interactions and they are forced to take on a value of pv/'2 at zq. In contrast 
to the isoscalar, the isovector densities exhibit much variation across the interactions and, 
for each interaction, the scaling needs to be exclusively due to an approximately invariant 
interdependence between the isovector and isoscalar densities and cannot be attributed to 
any imposed auxiliary condition. 

Figure M next compares the isoscalar and isovector densities in symmetric half-infinite 
matter, for the same sample interactions. For each interaction, the isovector density is 
generally pushed out relative to the isoscalar density. Amongst the interactions, for the Z 
interaction with the most negative L, though, the two densities are rather close to each other. 
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FIG. 6: Net density profiles at different asymmetries in half- infinite nuclear matter within SHFM, 
for sample interactions. For the SkI5 interaction, the asymmetry of 0.3 is already past the drip 
value of rjy = 0.288 and the corresponding profile is not shown. The r] = SkMP density is the 
same as the po-density in Fig. \5\ 

As L increases, the densities separate more and more, with the most pronounced differences 
between these densities appearing for the SkI5 interaction in the figure. The growing differ- 
ence with growing L again conforms with the expectations developed in Subsection IIIIDI 
Given the qualitative utility of the local approximation, it can be of interest to test that 
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FIG. 7: Profiles of the isovector density, Eq. ()44p . at different asymmetries in lialf- infinite nuclear 
matter within SHFM, for sample interactions. For the SkI5 interaction, the asymmetry of 0.3 is 
already past the drip value of rjy = 0.288 and the corresponding profile is not shown. The abscissa 
scale is numerically identical to that in Fig. [6j 

approximation quantitatively. For that purpose, in addition to the already displayed 77 = 
isoscalar and isovector densities, we plot in Fig. [8] the local approximation to the isovector 
density, of Eq. (1761) . It can be seen that the actual isovector density p°(-2) oscillates around 
the local approximation at the higher isoscalar densities, down to the net density reaching 
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FIG. 8: Comparison of isoscalar and isovector densities, and of the local approximation to the 
isovector density, in symmetric half-infinite nuclear matter within SHFM, for sample interactions. 
The longer- and shorter-dashed vertical lines for the specific interaction indicate, respectively, 
the location where the net density is equal to the quarter of normal density and the location of 
a classical return point for the Fermi wavevector directed along the z-axis, when the spin-orbit 
potential is disregarded. 
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about a quarter of the normal density, po/4, for the interactions in the figure. At fower 
net densities in the matter, the actual density and the approximation separate from each 
other, with the local approximation thereafter strongly overestimating p^ai'^)- follows 
that, for arguments r of the SHFM-operator S such that p(r) < po/4, this operator becomes 
significantly nonlocal. 

Increasingly nonlocal character of the operator iS, with a drop of density in the surface, 
is expected on the basis of the result fll04p . In that result, see fllOSp . the gradient corrections 
proportional to ti and ^2 represent effects of the range of interaction, while the corrections 
proportional to represent quantal nonlocalities. It is seen in the expressions (11051) that 
the fi^ correction-terms will dominate over the finite-range terms at sufficiently low densities 
and eventually also overwhelm the zeroth-order term. Overall, it appears that the nonlo- 
calities enhance the operator S in the far-out surface region making its inverse reduced and 
producing physically expected exponential fall-off for the isovector density pa- 

The close proximity of the isovector density pa to the local approximation, for the wide 
range of higher densities, opens up a chance of determining details in the density depen- 
dence of symmetry energy, from the systematics of proton densities alone. Correspondingly, 
it becomes of interest to understand better both the agreement of isovector density with the 
local expectation and the deviations. For this, we turn, in the following, to the Wentzel- 
Kramers-Brillouin- Jeffreys (WKBJ) approximation for the single-particle wavef unctions. 

E. WKBJ Analysis 

The wavelength for the oscillations of isovector density in Fig. [HI around the local ap- 



proximation, can be easily recognized as that expected for the Friedel oscillations |89|, 
\ = Ti /kp — 7r/(1.33 fm~^) = 2.36 fm. Corresponding oscillations can also be detected 
in the isoscalar density, but they are there of a considerably lesser amplitude than in the 
isovector density. Friedel oscillations generally arise when a disturbance of the system, 
forcing nonuniformity, such as the surface, acts to synchronize the density oscillations for 
individual component wavef unctions. Other than generating the oscillations, the specific 
synchronization of components due to the surface acts also to slow down the rise of density 
with distance away from the surface, when moving into the interior of matter. The Friedel 
oscillations are much stronger for the isovector density, due to the dominant contribution of 
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the states near the Fermi surface, to the difference of neutron and proton densities. 

To understand the oscillations mathematically, we invoke approximate WKBJ solutions 
to the SHFM equations fillip . In the classically allowed region, where the equation 

e,, = '-'^^^^^ ^ + U,{z) + X W,{z) , (141) 

has real solutions k^{z), the WKBJ wavefunctions are of the form 



TT 



^k^krxgi^) ^ sin i^j dz' kl^,^,^{z') + -j . (142) 

In the above, Zkj^k^xq is the classical return point ending the classically allowed region that is 
assumed to extend here up to infinity in the negative ^-direction, and Vz is the 2;-component 
of the velocity, Vz = h~^de/dkz = hkz/m*. In the classically forbidden region, where (114ip 
has purely imaginary solutions, the WKBJ wavefunctions are of the form 



,oo 



iPk^k^x,{z) ~ J^-^exp I - y dz' \kl^krxgi''')\ j ■ (143) 

Those wavefunction forms actually apply some distance away from the return point. Deep 
into the forbidden region, the wavefunctions are exponentially depressed and, in the con- 
siderations below, we will just put them to zero beyond the return point. With this, the 
density of nucleons q, at position z, becomes 

4 



4 



(27r) 



(27r)3 ^ 



/ dkz dk± sin^ (■ ■ ■ ) (^144) 
J dkz dkj 



1 + sin f 2 / dz' kl^^^^^{z' 



In obtaining the second to the last expression, we have converted the wavevector integration 
from one over the asymptotic to one over the local wavevector-components. It should be 
noted that those wavevector integrations are confined to a half of the Fermi sphere. In ar- 
riving at the last expression, we have used the identity sin^ (a + ^) = ^ + ^ sin (2a). Effects 
of the spin-orbit term on nucleon densities cancel out to the lowest order. For that reason 
and also to make the local Fermi sphere isotropic, those effects will be further disregarded. 
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l + sin 2 / dz' kl^^^^iz') 



In consequence, the density of nucleons q becomes 

With the disregard of the spin-orbit term, the local Fermi wavevectors in (11451) follow from 

If we were to consider the general case of mass potential B, we would need to decompose 
the local kinetic energy Tg in a similar manner to f ll45p . However, since our primary goal is 
to gain an insight, we shall confine ourselves to the case of m* = m, avoiding the need for 
specific decomposition. 

Mathematically, the Friedel oscillations can be tied to the sine term on the r.h.s. of fll45p 
or fll44p . Integration over wavevector components for that term leads to an averaging of 
the sine over argument values. While, in general, wide-range variations of the argument will 
make that average approach zero, any residue of averaging is likely to reflect possible nonuni- 
formities in the integration, as far as values of the sine argument are concerned. An obvious 
discontinuity is the termination of integration at kpq, producing oscillations of the form 
sin {2kpq z + 6). Otherwise, the local amplitude of oscillations depends on the distance from 
the surface. The larger the distance, the faster will be the variation of the argument of sine 
in f ll44p . with the variation of wavevector, and the more suppressed will be the contribution 
of the sine to a density. Furthermore, different amplitudes will emerge for different types of 
densities. Thus, isoscalar density involves 3-dimensional integrals over parameters affecting 
the sine argument, cf. (11450 . However, isovector density involves differentiation of nucleonic 
densities with respect to asymmetry, cf. fH6|) . which produces a sine-function integral reduced 
in the dimension to 2. In effect of the reduced averaging of the sine function, the Friedel 
oscillations turn out to be significantly larger, and they fall off slower with distance, in the 
isovector than in the isoscalar density. 

Considering the vicinity of zero asymmetry for the functions, we now represent nucleon 
densities as 

P,W = %^(l + ^.(^)), (147) 
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where 

r. I I Ai. 10/ A J I 



2ti k'^pq{z) Jk<kFq{z) 
The net nucleoli density at 77 = is then 



Qii^) ^ o_ J /.A / dfc Sin ( 2 / dz' kl , ,^(z') ) . (148) 



2k%{z) 



PoW = ^fV^(l + ^W). (149) 
We write, moreover, the /ia-derivative of pnp at = 0, for use in the isovector density (jl6l) . 



as 



where, with fll47p . 



On evaluating the /iQ-derivative of both sides of Eq. fll46p for neutrons, and on using 
Eqs. (UnD, (EHD, (El]), ffTOTH . ffT49l) and ffTHni) . we arrive at the following expression for the 
isovector density 

pl-iA .l., «^(l + ^^V (152) 



S{p)^^!^^Q^ Sip) \ QmS 

The last approximation is valid at large distances from the surface, in the classically allowed 
region, where we expect 1 ^ T ^ Q . 

Equation (I152p represents the isovector density in terms of the local approximation to 
that density combined with an oscillatory correction. Coefficient for the correction, whether 
examined in the denominator of the middle expression or on the r.h.s. of fll52p . involves the 
ratio of the kinetic contribution to symmetry energy to the net symmetry energy. Corre- 
spondingly, the Friedel oscillations are expected to be stronger for interactions with lower 
symmetry-energy values around the normal density, which indeed appears to be the case for 
the examples in Fig. [HI For the simplistic model of a step-like barrier at Zq, representing 
a rapidly rising isoscalar optical potential f/, the functions Q and may be arrived at in 
the following analytic forms: 

3 



Akl (zo - zY 



sin(2fc^ {zq - z)] 
cos{2k,izo-z)) . _ . 



(153) 



^(.)^-^^-^ sin (2^.(^0-.)). (154) 

Following those results, oscillations in the isovector density, governed by JF, should die out 
rather slowly with the distance {zq — z), as distance inverse, which appears to be borne out 
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by the results in Fig. [HI On the other hand, oscillations in the isoscalar density should die 
out rather quickly, as distance inverse squared, consistently with the previous qualitative 
expectations and with practical findings. 

A couple of other observations can be made in connection with the considerations above. 
Thus, it is apparent that the region of potential validity of the local approximation in the 
1] ^ matter is necessarily limited by the farthest lying classical return point for symmetric 
matter. At the general level, this criterion is consistent with the expectation that the range 
effects must be ignorable for the approximation to apply. An interesting aspect of this 
criterion is, however, that it is independent of the symmetry energy and tied only to the 
features of symmetric matter. Indeed, in spite of the widespread variation of the symmetry 
energy in Fig. [HI the position where local approximation breaks down, in terms of net density, 
remains approximately the same. The farthest return point is of course that for a Fermi 
wavevector directed along the 2;-axis and, for reference, we indicate the location of those 
points in Fig. [8l The return points for lower wavenumbers or for wavevectors not pointing 
along the z-axis are closer to the interior of the matter; the coarse end-location for the 
validity of the local approximation, of isoscalar density being about po/4, represents the 
return point for a Fermi wavevector directed at about 35° to the 2;-axis. 

Another observation, related to the preceding considerations, pertains to the rj- 
systematics of the tails of isoscalar and isovector densities. The relative variation of density 
in those tails appears to be significant with higher t] in Figs. [6] and U\ for at least two of the 
sample interactions, SkMP and SkI5. Features of the tails combine the effects of symmetry 
energy close to the classical return points and the effects of chemical potentials on the for- 
bidden region. Figure [H which displays details of different densities on a logarithmic scale, 
at different rj for two of the sample interactions, MSk9 and SkMP, can serve as an illustration 
in the discussion of understanding of the tails. At p ^ Po/4, larger deviations of symmetry 
energy from , for SkMP than for MSk9 interaction, and larger dS/dp, yield both stronger 
changes in the isoscalar p(r) with r]"^ for SkMP, as expected from Eqs. (1391) and ( |78i) . Changes 
within the second order in one of the densities, such as for SkMP, generally impact the other 
density within that order. Regarding the forbidden region, insights come from the behavior 
of the WKBJ wavefunctions far into the forbidden region, described in Eq. fll43p . With 
integrations in the wavevector space over squared wavefunctions producing densities, it is 
apparent that the neutron and proton densities need to fall off primarily exponentially in 
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FIG. 9: Densities on a logarithmic scale, in the surface region of half-infinite nuclear matter, for 
sample MSk9 and SkMP interactions in the top and bottom panels, respectively. The abscissa 
scale is numerically identical to that in Figs. El [7] and El The left panels show the approximately 
invariant densities, isoscalar (solid lines) and isovector (dashed), with the order from left to right 
in the tails corresponding to r/ = 0, 0.1, 0.2 and 0.3, respectively. The right panels show proton 
densities (dashed lines) for, respectively, rj = 0.3, 0.2 and 0.1, from left to right in the tail, and 
neutron densities (solid) for, respectively, rj = 0, 0.1, 0.2 and 0.3, from left to right in the density 
tail. 
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the far-out forbidden region, as exp(—2Kpq z) = exp{—z/dq), where upg = ^ —2m Hq/h and, 
thus, dq ~ 0.57 fm for ^q ~ 16 MeV. For > then, on account of the neutron contribution, 
with neutron chemical potential being greater or equal to the proton potential, the isoscalar 
and isovector densities should fall off as exp(— 2-^/— 2m finz/h). Given a milder neutron than 
proton slope for rj > 0, in fact, sufficiently far out, the neutron density needs to eventually 
strongly dominate over the proton density. Looking at the tails of proton and neutron den- 
sities in Fig. [9l far from the surface, it is seen that these densities indeed eventually begin 
to fall off exponentially, by about e per 0.6 fm at ?7 = 0. As has been anticipated, the fall-off 
slopes are always about the same for p and pa at a given rj, and, ai rj > 0, about the same for 
these two densities and p„. As neutron chemical potential increases with an increase in rj, the 
slope of fall-off of the densities becomes milder. When one moves then farther and farther 
out of the matter, for ?7 > 0, the accumulated effects of changed density slope, while possibly 
small on absolute scale, can become arbitrarily large on the scale of density scaled with the 
density for rj = 0. However, besides slopes, the behavior of densities in the allowed region 
matters as well for the forbidden region, in setting of the initial conditions for the fall-off. 
Incidentally, as critical asymmetries nearly coincide for MSk9 and SkMP, see Table [H there 
is a similarity in evolution of pn with changes in rj, for the two interactions. That similarity 
implies a degree of similarity in the slopes alone, that can be checked in Fig. |9l 



F. Coefficients Characterizing Nuclear Surface Energy, as and af 

For a self- sustained system in the absence of long-range interactions, the energy associated 
with a boundary can be simply defined as the difference between the actual energy and the 
energy expected if the boundary were not there. 



E, 



A 



dr 



(155) 



The last expression has been arranged to make the integration over depth within the matter 
finite in the limit of A — > oo, when the net energy can be represented in Skyrme form. In 
the case of half-infinite matter, Eq. f ll55p yields surface energy per unit area of the form 



^Skyrme (^) 



(156) 



From the energy per unit area, principally, the coefficients as and af in an energy formula 
may be inferred. Specifically, with the area of S ~ 47r Tq A^/^ for a large spherical nucleus 
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at — > 0, we find, from (fTTi) . tliat tlie coefficients are related to the energy of lialf- infinite 
matter with 

(157) 

In practice, some problems emerge when trying to utilize Eqs. (11560 and fll57p directly for 
determination of the coefficients characterizing nuclear surface. 

Thus, while the subintegral function in fll56p nominally approaches zero with increasing 
depth into the matter, i.e. z — oo in our convention, in practice it is not possible to em- 
ploy (11561) down to any depth in matter. This is because errors on the values of e{z) and p{z), 
from numerical integration of the SHF equations, do not vanish with increasing depth. If the 
use of (I156P were insisted upon down to any depth in the matter, the accumulated effects of 
errors across the volume would eventually overpower any surface contribution to the sought 
integral. Outside of the matter, errors from the integration of the SHF equations surge, 
at least relatively, close to the start of integration. In consequence, in integrations across 
surface of half-infinite matter, when seeking or other surface quantities, we employ 

a profile function P{z) to control errors accumulating with distance away from the surface, 
by choosing to approximate integration such as in (11560 with 



§ ~ / dzP{z) 



The profile function takes on the values 



eskyrmc(2;) " ( ^ ) P(^) 



V 



;i58) 





at 






at 


Z - Zq\ 



P{z) = { . (159) 



where zq represents surface location. We further make P{z) vary smoothly between zp and 
zp + Azp aiming at a reduction of the error in integration associated with an interplay 
between the termination of integration and the Friedel oscillations in nucleon and energy 
densities inside the matter. We typically employ zp ~ 6fm and Azp ~ 6fm. Results that 
will be presented, obtained using P, are stable with respect to variations in zp and Azp 
around our choices. 

When trying to extract the af coefficient from Eq. (11581) . an additional to the above 
problem arises, particularly serious for interactions with negative L. Thus, when extract- 
ing af , the energy values at small rj must be considered. This makes the contribution of 
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surface-symmetry energy to the net energy to be of net third-order in smallness, with one 
order due to the surface-to- volume ratio and with the two remaining orders due to the asym- 
metry. For a negative L, the surface contribution to the net symmetry energy is particularly 
small. The problem can be circumvented by employing Eq. flHTl) for calculating af , rather 
than (11561) . Namely, upon adding and subtracting A/a^ from the right-hand side of (1511) . 
for large but finite spherical system of mass number A, we find 

From the above, we find that the surface symmetry coefficient may be expressed as 



1 47rr, 



af a 



-^2 r 

^ dz {pUz) - p'{z)) , (161) 

a J 



where the r.h.s. integration is over the direction perpendicular to the surface. In practical 
calculations with Eq. (11611) . we also employ a profile function such as in (11581) . Reaching 
a desired accuracy for af with ffT6TD requires operating within two orders of smallness, 
rather than three, due to surface-to- volume ratio and due to difference p„ — pp, of the order 
of 7], in calculating p°. With the typical values of asymmetry that can be afforded in our 
extractions of af, r] ^ (0.03-0.15), the error on turns out actually to be (1-2) orders 
of magnitude smaller when using (I16ip than when using (I156p . In the following, we will 
extract the coefficients as and af , relying on the formulas above, and we will compare t 
results to those of preceding Hartree-Fock calculations for half-infinite matter by Kohler 

7, 8, 9^. Also, we shall compare the extracted coefficients to 



ne 
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961 for large, up to A ~ (10^ - 10^), 



and by Farine, Pearson et al. 

those deduced from spherical Hartree-Fock calculations 
systems. 

Interestingly, our Eq. 01611) is equivalent to a droplet-model expression considered in the 

n n n 

context of their calculations by Farine, Cote and Pearson [7|, |8| (see |3[ for a functional 
derivation by Farine), and ultimately exploited 9l(| in the context of Farine's code, where 



«f 2a^ ro vv^o r]v J { pl^ p^ 

The equivalence of the results (I162p and (I16ip can be seen by noting that at low volume 
asymmetries, rjy 0, the asymptotic nucleon densities become p^ ~ 3(1 ± ?7v/2)/(87rrQ), 
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FIG. 10: Different energy densities in the surface region of symmetric lialf-infinite nuclear-matter, 
for sample Skyrme interactions. Shown are the Skyrme energy density in the standard SHFM 
calculations, egj^yj-j^e Eq. (j80|) . the energy density from the calculations with spin-orbit term 
ignored, ^skyime' energy density in the Thomas- Fermi approximation with gradient terms ig- 
nored, Cgj^y^jj^g, and, finally, the energy density expected if the energy per nucleon were constant, 
{E/A)v po- The abscissa scale is numerically identical to that in Figs. El [7] and [51 



allowing to rewrite the r.h.s. of fll62p as 



V 



a 



i^il lim 



dz 



Pn{z) pp{z) 



dz 



l + rjv l-rjy 

Pn{z) - pp{z) 



Vv 



(163) 



where the r.h.s. yields next the r.h.s. of fll6ip . cf. Eq. (H6l) . Notably, because of the Friedel 
oscillations, Farine et al. ended up never employing Eq. (11621) directly in their analyses. 
Instead, they have examined the separation between neutron and proton surfaces at finite rjy, 
with which separation they have replaced the integral that appears on the r.h.s. of Eq. 
With this, they proceeded along the lines of the droplet model 27 1, see also 



Figure [To] displays different energy densities in the surface region of symmetric half- infinite 
nuclear matter, for two sample Skyrme interactions, MSK9 and SkMP. Those densities 
include the two energy densities, eg]^yj.jjjg and {E/A)y pq, of which the integrated difference 
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yields the surface energy fll56p . and they further include the energy density Cgjl^^'^g calculated 
with the spin-orbit term suppressed and, finally, the energy density calculated in the Thomas- 
Fermi approximation, Cgj^^j^^. In obtaining the last energy density, the density of kinetic 
energy is taken in the local approximation of Eq. (jHS]) and the gradient terms in any density 
are ignored. 

For moderately subnormal nucleon densities in Fig. [TOl the energy densities {E IA)v Po{z) 
and e^^yj.^^{z) are not much different from each other. This just reflects the fact that the 
energy per nucleon minimizes in uniform matter at normal density and, in consequence, it 
changes little with moderate deviations of density from normal. Across nucleon densities, 
the energy density e^^^yT-mei^) oscillates in the figure in the vicinity of e^^yj.^^{z). Dominant 
and positive contribution to the difference of energy densities e'l^y^^^ — {E/A)vPo, in the 
expression (11561) for the surface energy, is seen to come from the far-out tail of the density 
distribution, from essentially the classically forbidden region. In fact, in the outmost tail, the 
potential energy approaches zero faster than does the kinetic energy, making the net energy 
density even positive. The positive contribution to the difference in (11561) is next, in the 
matter direction, much reduced by the negative contribution from just moderately subnormal 
nucleon densities in the density tail. When moving further in position into interior of the 
matter, the energy densities in the difference in (11561) and, correspondingly, the difference 
itself oscillate in a Friedel fashion, with a wavelength of ~ n/kp ^ 2.4 fm. The precise 
location of maxima and minima in eskyrme, and in the energy difference, depends on the 
choice that was made for the kinetic energy density. Thus, when expressing the kinetic 
energy density in terms of a square of wavefunction gradient in (1821) . the kinetic energy 
oscillates approximately out of phase with respect to the density. If the kinetic energy were 
expressed, though, in terms of a laplacian, rather than gradient squared, the oscillations of 
kinetic-energy density would have been approximately in phase with the nucleon density. 

Surface energy coefficients, calculated by integrating the rj = difference of energy den- 



sities in (I158p . to obtain the energy for elementary surface area, 

as = ATTrl^, (164) 

are provided for various Skyrme interactions in Table [H The value distribution, across 
interactions, is rather narrow, with the average coefficient value of 17.6 MeV combined with 
the standard deviation of 1.1 MeV. 
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An interesting issue is how the surface coefficient may be affected by the effects of a spin- 
orbit couphng acting in the surface region. Looking at the Skyrme energy functional (l80l) . 
where nucleon densities multiply spin-density divergencies in the coupling, we see that this 
coupling allows a nuclear system to lower its energy by developing spin densities. Divergence 
of a spin density obviously integrates to zero. However, energy of a system can change when 
positive and negative divergence values are correlated with different magnitudes of nucleon 
densities. For positive 64-coefficient values in the functional flHOj) . the energy gets lowered 
when a spin density directed outward develops in the surface region. In that case, the di- 
vergence is positive in the inner part of the surface and negative in the outer allowing the 
energy density to drop more in the inner than in the outer part of the surface. The polariza- 
tion in the system is moderated, in particular, by the kinetic energy density that grows, but 
only quadratically, with spin density. For reference, in addition to other energy densities in 
Fig. [lOl we show there also the Skyrme energy density e'^^^^ obtained when suppressing the 
spin-orbit coupling. It is seen that, indeed, inclusion of the coupling lowers energy density 
more in the inner part of the surface than in the outer. Lowering of the system energy leads 
to lowering of the surface coefficient. Notably, inclusion of the coupling affects also {E/A)v p 
in the difference of energy densities for the surface energy (11581) . but not as strongly as eskyrme 
there, cf. Fig. [51 Quantitatively, for the surface coefficients of the interactions illustrated 
in Fig. [101 we find af=° = 18.7 MeV for both the MSk9 and SkMP, when the spin-orbit 
coupling is suppressed. With the coupling, we find as = 17.1 and 16.6 MeV, respectively, 
for the the MSk9 and SkMP interactions, see Table [H A ~ 10% magnitude in reduction of 
as is typical for the Skyrme parameterizations, when switching on the spin-orbit coupling, 
see also the forthcoming Table [ITl To our knowledge, Stocker ||88|] was the ffist to assess the 
sign and magnitude of the change in the surface symmetry coefficient when including the 
spin-orbit coupling in the surface description. 

Besides the values of as, Table [J gives, further, the surface symmetry coefficients af ob- 
tained by integrating the 1] = difference of isovector and isoscalar densities in Eq. (11611) . 
As might be expected from Fig. [HI where the relation between isoscalar and isovector den- 
sities varies widely from one interaction to another, the values of af vary widely between 
individual interactions, from 10 to 60 MeV. Figure [TT1 illustrates next the changes in surface 
energy (11581) of an elementary area 47rrQ, with changing asymmetry, for the sample Skyrme 
interactions. Magnitude of the changes in energy should be compared to the magnitude of 
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FIG. 11: Symbols represent change in the surface energy of an elementary area, from (jl58p . for 
different Skyrme parameterizations, plotted vs square of asymmetry: volume asymmetry in the 
left panel and effective volume asymmetry, ijy = /Xa/(2o^), in the right panel. Lines represent 
expectations based on Eq. (jl57p . with surface symmetry coefficients taken from Eq. (jl6ip . 



the surface coefficient fll64p . of the order of 17MeV for the Skyrme interactions. Note that 
the rise of surface energy with asymmetry in a finite nucleus would be accompanied by a drop 
in the volume symmetry energy, compared to a system with no surface, due to asymmetry 
moving out to the surface region. As evident in Fig. [HI the surface energy generally changes 
quadratically with asymmetry. At small asymmetries, barely visible deviations from the 
lines, representing expectations from Eqs. fll57p and fll6ip . stem from numerical inaccura- 
cies, primarily arising from the use of Eq. fllSSp to calculate the surface energy. At higher 
asymmetries generally true fourth-order terms in asymmetry come into play. In particular, 
some differences between the regularities in the left and right panels become apparent at 
higher rj and L; these are due to the differences that, depending on interaction, develop 
between rjv and riy at higher asymmetries, as a consequence of fourth-order terms in the 
energy of uniform matter. 

For the Skyrme parameterizations with 64 = 64, cf. Eq. (!84l) . the spin-orbit coupling 
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has even less effect on af than on as- Unhide as, the coefficient af increases rather than 
decreases with the switching on of the spin-orbit couphng, usually by less than 10% and 
often by much less. Notably, though, at constant volume asymmetry, the increase in af 
implies dropping of the surface energy, cf. fll57p . just as a decrease in does. Without the 
coupling, we find e.g. af^^° = 22.5 and 11.1 MeV, respectively, for the MSk9 and SkMP 
parameterizations, and af = 23.0 and 12.1 MeV, respectively, for these parameterizations, 
with the coupling. Additional exemplary results from 64 = 64 Skyrme parameterizations can 
be found in the forthcoming Table HTl For parameterizations with different 64 and 64, the spin- 
orbit coupling may have a stronger effect on the af coefficient, than for 64 = 64, seemingly 
due to a strengthening of the spin-isospin correlations in the surface region. E.g. for the 
SkO parameterization, we find = 13.2 MeV and af = 14.9 MeV, without and with 

the coupling, respectively. For SkI4, we find af = 17.5 MeV and af = 20.8 MeV. 



G. Comparison to Coefficient Values in the Literature 



Our results for the coefficients of surface energy are next compared in Table [TTl to the 
results obtained by others in the literature within the SHFM calculations of asymmetric 
systems. Specifically, our coefficients are compared to those obtained within the calculations 
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of semi-infinite matter, done by Kohler [5|] and by Pearson et al. 
94l |. The latter results stem predominantly from the code by Farine [95|. The dependence of 
the surface energy on asymmetry has been expressed in the past calculations of semi-infinite 
matter in terms of the coefficient Q from the droplet model, cf. Subsection III CI That 
coefficient differs from af by simple factor, cf. (fT6|) . allowing for a simple transcription of 
the results obtained for semi-infinite matter. The calculations by Kohler were done without 
spin-orbit coupling and, in comparisons to his results, we have suppressed that coupling 
as well. In addition, in Table [ITl our coefficients for the surface are compared to those 
deduced in the calculations of spherical nuclei with large mass, up to A ~ (10^ — 10^), by 
Reinhard et al. [96]. Those authors have assessed the coefficients of expansion of energy in 
powers of y4~^/^ and rf, for different interactions. If we expand the symrnetry coefficient ([1] 
in A~^/^, we find that we can compute af from the coefficient flssym 

... «f 



m 



96], with 



(165) 



"ssym 
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As will be mentioned later, in addressing the likely coefficient values, the expansion of the 
symmetry coefficient f|T^ lacks justification for nuclear masses encountered in nature. 



TABLE II: Comparison of coefficients characterizing nuclear surface energy, 
from SHFM calculations of nuclear systems by different authors, with sources 
for outside results provided in the superscripts. The surface symmetry coef- 
ficient of is related to the coefficient Q, stemming from the droplet model, 
with of = 4Q/9, and to coefficient Ossym iu the large-^ expansion of nuclear 
energy [OQ], with of = -{a^f /a^syra- 







Value (MeV) 


Name 


a 


Kohler 


Pearson et al. 


Reinhard et al. 


Our 


SII 


51^=0 
"a 

as 
aa 


19.8^5] 

ii.et^ 






20.20 
17.27 
19.38 
17.54 


SIII 


"a 

as 

"a 


19.6^5] 


19.9l^ 
23.1^8] 




19.85 
21.79 
18.54 
21.77 


SIV 


as 
a'a 




19.6l^ 

15.1^81 




18.87 
13.39 


sv 


„sw=o 

"a 

as 
aa 


19.6^^1 


19.8^8] 
11.1^81 




20.44 
9.95 
19.16 
10.27 


SVI 


as 

"a 




18.3ia 
26.7ta 




18.10 

26.27 


SkMs 


as 




18.0^90] 




17.46 




"a 




26.7f9Q.] 


174f96] 


14.48 



Table Continued on Next Page. 
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TABLE II - Continued 



Name 


a 


Value (MeV) 


Kohler 


Pearson et al. 


Reinhard et al. 


Our 


SKa 


„sw=o 

Q, 


19.8^^1 
9.3^ 






19.97 
12.06 


SKb 


„sw=o 

CL 


19.8l^ 

7.6^5] 






19.97 
10.03 


SkP 


as 
af 






18.2t^ 

20.0^96] 


18.18 
18.02 


SkI3 


as 
al 






18.0126J 
16.2[26J 


17.77 
12.77 


SkI4 


as 
al 






177[96] 
25.6t26J 


17.48 
20.83 


SLy4 


as 
a'a 






18.4t26J 
19.0l26j 


18.24 
16.60 


SLy6 


as 
af 




17.74f9iJ 
14.8f9iJ 


17.7^9^1 
20.0^96] 


17.53 
17.12 


SkO 


as 
af 






17.3t^ 
17.6[96] 


17.14 
14.94 


BSkl 


as 
an 




17.54i2^ 
20.3i2^ 


17.5^96] 

2i.5EeJ 


17.22 
23.72 


BSk2 


as 
af 




17.54[7i 91, 92] 

20.4^911, 30.2^71' 




17.14 
23.15 


BSk3 


as 
af 




17.5[71,92] 

29.8Ei'^ 




17.16 
23.37 


BSk4 


as 
al 




17.3il^ 
33.8E1] 




16.91 
23.01 


BSk5 


as 
al 




17.5il^ 
23.li2^ 




17.04 

20.65 



Table Continued on Next Page. . . 
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TABLE II - Continued 







Value (MeV) 


Name 


a 


Kohler 


Pearson et al. 


Reinhard a/. 


Our 


BSk6 


as 






173[96] 


16.74 




«f 




I9.<m, 23.7il^, 36.9^^ 


23.8t26J 


22.73 


BSk7 


05 




17.3Ei,9^ 




16.70 








20.1^91], 35.6i2^ 




22.81 


BSk8 


as 




17.64t2^'-2i'^ 




17.12 




af 




20.2^'^, 24.0iZ2] 




23.13 


BSk9 


05 




17.92t2i'93] 




17.42 




«f 




15.8[9i,93] 




17.83 


BSklO 


o-s 








17.51 








15.6^73] 




17.48 


BSkll 


o-s 








17.32 




af 




15.6iZ2] 




17.42 


BSkl2 


05 




17.7f73] 




17.30 








16.0^73] 




17.51 


BSkl3 


05 








17.34 




a^ 




15.6Ea 




17.34 


BSkl4 


as 




17.6^' -2^ 




17.17 




4 








17.22 



When carrying a case-by-case examination of the entries in Table [IT| it becomes apparent 
that results for the surface coefficient as agree fairly well between different authors. As to 
results for the surface symmetry coefficient af , a degree of agreement is found in many cases 
but in a number of other cases there is a disagreement by factors even in excess of 1.5. 
The status of the two coefficients in the literature is further illustrated with two panels 
in Fig. [121 For different Skyrme parameterizations with or without spin-orbit coupling, 
the values obtained in the literature are plotted in those panels against our values. In case 
of an ideal agreement, the results should line up with the diagonal lines in the panels. 

68 



In assessing the level of agreement between the results in literature, on the basis of the 
figure, it should be noted that the range of values displayed in the panel for as is lower by 
one order of magnitude than in the panel for af . 

Regarding the three author groups, we appear to have the best overall agreement with 
Reinhard et al. [96]. Their results, though, are also later than most of other results that we 
compare against. On the average, their results for as are higher by just ~ 0.25 MeV than 
ours for the coefficient. The results by Pearson et al., for the surface coefficient, appear 
systematically shifted by a bit more, ~ 0.4 MeV, which still represents just 2% of coefficient 
value. By contrast, in the overall assessment, the agreement between the results for af 
turns out to be somewhat unsatisfactory, given the discrepancies encountered that are of 
the same order as coefficient values. In practice, this demonstrates the difficulty in working 
with contributions to system properties that are small due to their simultaneous association 
with surface and with asymmetry. That difficulty is further underscored by the fact that 
Pearson et al. have quoted significantly different values characterizing surface symmetry 
energy, for the same interaction, in different publication years, cf. Table [TTl 

Realizing the fragility of their results for the coefficient characterizing surface symmetry- 
energy, Pearson et al. have eventually reassessed the procedures for extracting the coefficient 
values, specifically in Appendix A of Ref. [91|]. They have examined there the stability of 
extracted coefficient values, when changing asymmetry in the extraction. From the three 
tested procedures tested, they have found the one based on Eq. fll62p to be the most stable 
and, thus, preferred. Note that before 7|, for the sake of result stability, Pearson et al. have 
found it beneficial to approximate the integral on the r.h.s. of fll62p . with a displacement 
of densities. While there is a significant improvement in the a gree ment of our results with 
those of Pearson et al. following the reassessment 0, [91I, 91, 94], disa gree ments between 
the two sets still can reach the order of 20%, cf. Table [TTl Reinhard et al. [96] have employed 
a strategy close to that in Eq. fll57p and disagreements between their and our af -coefficient 
values reach ~ 25%. 

Our own tests of the stability of our results for af from Eq. f ll6ip . with respect to 
different technical aspects of calculations, indicate that errors of those coefficient values do 
not exceed 0.5%. The aspects of calculations that have been varied include the already 
mentioned mesh in wavevector space, step of integration in space, and the profile function in 
space. In the context of the investigations in 9l(], it may be still worthwhile to compare the 
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FIG. 12: Comparison of the values of coefficients characterizing nuclear surface 
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91 



92 



93, 



energy m our 



94j, 



gg]. The 



calculations and elsewhere in the literature 
first and second panels show, respectively, the comparison of surface coefficients as and of surface 
symmetry coefficients af. Respective coefficient values are given in Table Ull and the comparisons 
of coefficients with and without spin-orbit coupling are lumped together. In the case of an ideal 
agreement, the results, as plotted, should line up with the diagonals in the panels. 
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results on af obtained with different methods. We have mentioned that errors for a method 
involving dividing out two powers of rj, can exceed errors for a method involving dividing out 
just one power, by (1-2) orders of magnitude, for the small values of rj we employ. Partly, 
this can be compensated by employing a higher r] in extracting af with the inferior method: 
when one is resigned to a larger error, effects of any anharmonicity of the symmetry energy 
in 7] matter less. Overall, though, when the quality of calculations increases, errors decrease, 
whether for the more favored or the more inferior method of af extraction. Thus, differences 
between the results obtained with different methods can be exploited in general assessment 
of calculations. 

As to the intrinsic comparisons. Fig. [HI for one, can be used for testing visually the 
consistency between Eqs. f ll57p and (11611) in our calculations, for the several interactions 
in that figure. Otherwise, Table IIIII shows the values of af extracted with three different 
methods from our half- infinite matter calculations, when employing Eq. (I16ip . Eq. (I157p 
and also Eq. ( fTTl) . The last two methods both involve dividing out factors of the second 
order in asymmetry and, thus, are both expected to provide inferior results to the first 
method that is our standard. Our calculations for the table have been done for those 
Skyrme parameterizations for which Pearson et al. have also done calculations with different 
methods. Their results are also shown in Table IIIII for comparison. The largest deviation 
between the preferred and an inferior method in Table IIIII is 7% for our calculations and 
34% for Pearson et al. 

H. Further Discussion of the Isovector Density and of the Coefficient 

As we have shown, in the classically allowed region of the surface, p ^ Po/4, the isovector 
density pa follows closely the expectation of Eq. (1761) from uniform matter, with pa inversely 
proportional to the symmetry energy S{p) at a given net density p. The form of the adher- 
ence to the expectation, and the partial consequences of that adherence, are the following 
for the specific limits of S{p). When the symmetry energy remains significant at subnormal 
densities, for low positive or even negative values of the dimensionless parameter-ratio L/a^ , 
the isovector density pa stays close to the isoscalar density p over a significant portion of 
the classically allowed region of nuclear surface, as expected from Eq. ( 1761) . Outside of the 
classically allowed region, both densities drop rapidly to zero, in an approximately expo- 
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TABLE III: Values of surface symmetry coefficient, extracted following different indicated 
equations from the SHFM calculations of half-infinite nuclear matter, done for the different 
indicated Skyrme parameterizations. The results that were ultimately decided to be superior by 
the respective authors are represented in the roman font and other results are represented in italic. 



Name 


af [MeV] 


Our 

Eq. (dH) Eq. dlMJ) Eq. (jl57|) 


Pearson et al. 
Eq. (062])" Eq. dlSTj) Ref. 


SII 


16.63 17.54 18.26 


15.1 17.8 [8] 


SIII 


20.83 21.77 22.75 


19.1 23.1 [8] 


SIV 


12.88 13.39 U.27 


12.9 15.1 [8] 


sv 


10.40 10.27 10.80 


8.9 11.1 [8] 


SVI 


26.13 26.27 26.66 


32.4 26.7 [8] 


BSk8 


22.80 23.13 23.32 


20.0 13.3 [91j 



"All additional approximation is employed 7[. 



nential fashion. On the other hand, when, for high L/a^-ratio, the symmetry energy drops 
rapidly with density at subnormal densities, the isovector density pa remains significantly 
higher than p, across the classically allowed region, cf. Fig. [HI For the exponential fall-off 
of the densities in the classically forbidden region, the starting value for pa is then further 
significantly higher than for p, extending the region where pa dominates over p, cf. Fig. [91 
Note that, while L/a^ ratio is suitable for characterizing the shape of 5'(p), due to small 
variations of relative to variations of L between different interactions, the relative magni- 
tudes of L and the relative magnitudes oi L/a^ are usually interchangeable when comparing 
different interactions. 

In Table [Tj different parameters can be found, quantifying characteristics of the densities 
Po and p°, including respective diffuseness parameters, and (f^. The isoscalar diffuseness 
do is of obvious importance as quantifying the pace, independent of the symmetry energy, 
at which the isoscalar density approaches zero. We define each diffuseness parameter in 
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terms of the derivative of respective density at half of asymptotic value: 



Mo 



dpo 
dz 




dp! 



,0 



a 



(166) 



P0(2)=P0/2 



pO(2;)=Po/2 



We elect such a definition over the more common relation of diffuseness to the surface 
thickness as distance over which the density drops from 90% to 10% of asymptotic value. 
This is because of our intention to have consistent definitions for p and pa and because of 
the potentially complicated features of pa around normal density. The latter features are 
due to the relation of pa to symmetry energy and due to amplified Friedel oscillations in pa, 
as compared to p, cf. Fig. [HI For the isoscalar density, the definitions in terms of derivative 
and surface thickness yield nearly identical results. 

The isoscalar diffuseness parameters tend to be fairly consistent between different Skyrme 
interactions, averaging at 0.541 fm in Table [H with an rms deviation of 0.031 fm. When 
a larger deviation from the average is found for an interaction in Table [H a correlated change, 
compared to other interactions, is typically found to occur in other distance parameters for 
the surface. For the lowest L- values, for which the po/2-position of the isovector density lies 
within the classically allowed region, the isovector diffuseness needs obviously to depend on 
the shape of the dependence of symmetry energy on density, cf. Eq. fl76|) . At any fixed shape 
for the symmetry energy, such as characterized by the ratio L/a]^, though, the isovector 
diffuseness should further scale in proportion to the isoscalar diffuseness that sets the scale for 
variation of isoscalar density in space. Figure [13] shows the correlation between, on one hand, 
the isovector diffuseness from the SHFM calculations, scaled with the isoscalar diffuseness, 
and, on the other, the slope parameter L scaled with . A quite tight correlations is 
observed between the scaled parameters, with a linear fit to the results yielding 



Interestingly, the tight correlation extends well into the region of the isovector po/2-point 
lying beyond the po/4-point in isoscalar density, representing the start of a classically for- 
bidden region. Figure [T^ shows also the correlation between the isovector diffuseness and 
slope L, when the parameter scaling is removed. The lack of scaling results in a deterioration 
of the correlation; analogous deteriorations may be found for other tight correlations in the 
paper, when scalings implied by the physics are removed. 

More significant than the difference in diffusenesses, between the isovector and isoscalar 



< -do (1-062 - 0.117L/a^). 



(167) 
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FIG. 13: Correlation between the diffuseness of isovector density in symmetric half- infinite 
matter and the slope parameter L of symmetry energy in uniform matter. In the first panel, 
the isovector diffuseness df^ and the slope parameter L are scaled, respectively, with the diffuseness 
do of isoscalar density and with the value of symmetry energy at pQ. No scalings are applied 
in the second panel. Symbols in the panels represent results for Skyrme interactions in Table [H 
Solid lines represent linear fits to the results for Skyrme interactions. The dashed vertical line in 
the first panel separates coarsely those interactions, to the left of the line, for which the /3o/2-point 
in isovector density appears in the classically allowed region, from those interactions, to the right 
of the line, for which the /3o/2-point appears in the classically forbidden region. 

densities, is the overall displacement of those densities relative to each other. This displace- 
ment can be quantified in different ways. One possible quantification is in terms of the 
displacement A/2° of the po/2-density points for isovector and isoscalar densities, at r] = 0, 
see Table HI Another possible quantification is in terms of the ratio of symmetry coeffi- 
cients /a^, proportional to the integral over the difference of those densities at rj = 0, 
in the direction perpendicular to nuclear surface, cf. Eq. (11611) . To the extent that the 
densities could be both described in terms of Fermi functions or, otherwise, were simply 
translations one of another, the ratio of coefficients could, following (116ip . be approximated 
with 

(168) 



To 
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FIG. 14: Correlation between, on one hand, the relative displacement AR^ of isovector and isoscalar 
densities in half-infinite symmetric matter, from the /9o/2-points, and, on the other hand, the 
effective displacement Ag-R = (o^/af ) (ro/3). Symbols represent the variety of Skyrme interactions 
from Table [H The diagonal dashed line serves to guide the eye. 



Following fll68p . we define an effective displacement in terms of the coefficient ratio as 



3 a 



s 



(169) 



which is provided in Table [T] in addition to Ai?°. Surface displacement AR^ is found to 
range in the table from 0.21 to 1.13 fm. Figure [14] shows the correlation between the two 
displacements, for the Skyrme interactions. Over much of the range of variation, the dis- 
placements are rather close to each other. Only for the highest displacement values, Ae-R 
tends to prevail over Ai?°. 

Judging from the local approximation, and from Fig. [HI the most significant differences 
between the isovector and isoscalar densities should be associated with the highest L-values. 
Indeed, high L-values imply a rapid drop of symmetry energy at subnormal densities and, 
correspondingly, high isovector compared to isoscalar density values. Correlation between, 
on one hand, the symmetry parameter ratio a^/af , proportional to the integral over density 
difference, cf. Eq. (116 ip . and, on the other hand, the scaled slope parameter L/a^ , is shown 
in Fig. [151 The lowest values of the symmetry parameter ratio, of the order of 0.4, are 
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FIG. 15: Ratio of symmetry energy coefficients a^/of , for different Skyrme interactions, plotted 
vs the slope parameter L of the symmetry energy, scaled with the value of symmetry energy . 

indeed found for the lowest values of scaled L, and the highest values of the ratio, of the 
order of 3.4, are found for the highest values of L. In more detail, the isovector density 
is expected to stay elevated, following the symmetry energy in the local approximation of 
Eq. (176!) . down to the isoscalar classical- return density of ~ po/4. With this, a correlation 
can be expected between the displacement Ai?° of the isovector from isoscalar surfaces and 
the value of symmetry energy at po/4. Figure [TBI shows the correlation between Ai?° scaled 
with isoscalar diffuseness and the symmetry energy at po/4, scaled with a^. The correlation 
is the tightest from among those explored here for surface properties. Linear fit to the 
correlation, represented by a solid line in Fig. [161 produces 

Ai?° ~ do (2.71 - 3.42 S{po/4:)/al) . (170) 

The dashed line in Fig. [16] represents prediction based on the simplified geometric consid- 
eration illustrated in Fig. [TTJ where the densities are taken to vary linearly with position. 
That consideration yields 

A/?° c^do + dl (ar/5(po/4) - 2) , (171) 

and, for Fig. [161 we use (i° ~ (io (0.336 — 1.336 S'(po/4)/a^) from fitting the correlation 
of (i° with symmetry energy. While the formula from the simple consideration begins to 
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FIG. 16: Displacement of isovector- relative to isoscalar-density in symmetric half-infinite matter, 
plotted vs value of symmetry energy at po/4. The displacement and the energy value at po/4: are 
scaled, respectively, with the isoscalar diffuseness and the symmetry-energy value at pQ. Symbols 
represent results for Skyrme interactions from Table [H Solid line represents a linear fit to the 
Skyrme results. Dashed line represents predictions of formula (jl7ip based on simplified geometric 
considerations. 

overestimate the shift for the symmetry energies with a strongest density- dependence, it 
is apparent that the simple consideration grasps the essence of impact of the low-density 
symmetry-energy on the shift. 

Rounding up the set of correlations between surface characteristics and features of sym- 
metry energy in uniform matter. Fig. [TS] shows the displacement of the isovector- relative 
to isoscalar-surface for the Skyrme interactions, as a function of scaled slope parameter. 
Quality of this correlation does not change significantly when the shift is scaled with diffuse- 
ness, so we present the correlation without such a scaling, as possibly simpler to use. With 
the correspondence between the displacement AR^ and the effective displacement AgR pro- 
portional to /a^, cf. Fig. [HI Fig. [18] is a counterpart of Fig. [151 Linear fit to the Skyrme 
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FIG. 17: Geometric construction behind Eq. (jl71|) . In the surface region, the dependence of 
isoscalar density and of the isovector density on position, beyond the breakdown of local approxi- 
mation for is approximated in a Unear form, with the slope expressed in terms of the respective 
diffuseness. 




-2 2 4 



FIG. 18: Displacement of isovector- relative to isoscalar-surface in semi-infinite nuclear-matter, 
plotted vs normalized slope-parameter of symmetry energy. Symbols represent results for the 
Skyrme interactions of Table [H The line represents a linear fit to the Skyrme results. 
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results, which is indicated in Fig. [18], produces 

Ai?° ~ (0.458 + 0.165 L/a^) fm . 



(172) 



I. Some Context 

In discussing the isovector and isoscalar densities in the preceeding subsection, we have 
concentrated on the shift of these densities relative to each other and on their diffusenesses. 
To the extent to which and po niight be approximated in terms of Fermi functions, these 
characteristics could suffice for describing the average sizes of asymmetry skins in nuclei. 
The skin sizes have been, so far, the prime aspect of nuclear density distributions linked in 
the literature to the density dependence of symmetry energy. 

Specifically, in the context of asymmetry skins, let us consider neutron and proton rms 
radii for a nucleus of mass A. From fl52l) . we find for the mean squared proton radius: 

{r\={r') + ^{{r')a-{r')) , (173) 

where 

A 



(r') = T / dVr^pir), (174) 



and 



JdVr^p^ir) 



^ J dVr'p^.iv). (175) 



fdVpair) N 

With the N ^ Z interchange on the r.h.s. of fll73p yielding the neutron instead of the proton 
mean squared-radius, the difference between nucleon squared-radii may be represented as 

{r\'{r\ = ^^^^{{r\-{r')). (176) 

From the above, we finally get for the size of asymmetry skin 

/,2u/2 _ ,^2.1/2 ^ A{N-Z) (r^). - (r^) 

It turns than out that the size of asymmetry skin directly probes the difference in mean 
squared radii for the isovector and isoscalar densities. When a density can be approximated 
in the Fermi form, the mean squared radius can be expressed in terms of the po/2-position 
combined with diffuseness. Looking at Figs. [6] and [71 it is apparent that isoscalar densities 
are more likely to be well described in terms of a Fermi shape than isovector, particularly 
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for the most extreme symmetry energies, characterized by either very low or very high 
slope-parameters L. 

Outside of the scope of this paper, analyses of data point to a symmetry energy in uniform 
matter that is characterized by a slope in the upper range of the values possible for Skyrme 
interactions, especially for the slo pe scaled with a^. Thus, e.g. the analysis of excitation 
energies of isobaric analog states |97|, [OSj yields independent values of and af . While 
the volume symmetry coefficient from this type of analysis, ~ (31.5-33.5) MeV, comes 
out quite in the middle of values found for the Skyrme interactions, the surface symmetry 
coefficient, af ^ (9.5-12) MeV, comes out right at the lower end of values encountered for 
the Skyrme interactions. The coefficient ratio from that analysis is in the range o^/af — 
(2.8-3.3). That ratio produces the effective surface displacement in the range of Agi? = 
(ro/3) (a^/af) ^ (1.06-1.26) fm. Moreover, Figs. \\M and [T5] yield the respective ranges of 
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98| 



Ai?° ~ (0.85-1.05) fm and L/a^ ~ (2.4-3.4) or L ~ (78-111) MeV. The analysis 
is relatively model-independent, provided curvature effects play little role for heavier nuclei. 
If the latter were not the case, though, a bit softer symmetry energy would need to be 
deduced. 

The large displacement for the isovector surface, inferred above, would give rise to large 
predicted asymmetry skins, cf. Eq. (11771) . towards the top of those that appear possible 



within the Skyrme-Hartree-Fock descriptions 
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991] . Notably, still higher sizes of asymme- 



try skins can be encountered in the relativistic mean-field descriptions of nuclear systems, 
along with higher values of the L-parameter than encountered for the Skyrme parameter- 
izations 99|]. Aiming at constraints beyond those in 97|, 1981], on the symmetry energy in 
uniform matter, we intend to develop strategies for the constraints, further in this series, 
exploiting available experimental information on nucleonic densities, along the lines of devel- 
opments of the present paper. Otherwise, the isovector density, as a fundamental quantity, 
deserves attention on its own when analyzing data on nucleonic densities. 



V. CONCLUSIONS 



In this paper, we have considered nuclear energy in the macroscopic limit for a nucleus. 
Further, we have considered the Hohenberg-Kohn functional for a nuclear system, in terms of 
proton and neutron densities. Finally, we have carried out Skyrme-Hartree-Fock calculations 
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of half-infinite particle-stable nuclear-matter. In each case, we have concentrated on the role 
of neutron-proton asymmetry in the system and on the symmetry energy. 

In discussing nuclear energy in the macroscopic limit, we have shown that surface 
symmetry-energy emerges, as an unavoidable ingredient of the net nuclear energy, from 
simultaneous considerations of nuclear surface and symmetry energies. A consequence of 
this conclusion is that the net nuclear surface and volume energies combine in the similar 
manner as energies of connected capacitors in electrostatics. The role of charge flowing and 
distributing itself between the capacitors, is taken over by the asymmetry that distributes 
itself between the nuclear interior and surface, in proportion to volume and surfaces capac- 
itances. The capacitances for the interior and surface are proportional, respectively, to the 
volume and surface area and are inversely proportional to the volume and surface symmetry 
coefficients. 

When considering a continuous limit of the Hohenberg-Kohn functional, we have broken 
up the nuclear part of the functional into the functional for symmetric matter and a sym- 
metry term that is bilinear, in the lowest order, in the difference of neutron and proton 
densities. The kernel S in the symmetry term is generally nonlocal and depends on the 
net density. In the limit of weak nonuniformities, the kernel may be approximated in a 
local form as a (5-function multiplying symmetry energy in uniform matter, divided by lo- 
cal density. We have shown that, up to Coulomb corrections and terms of second order in 
asymmetry, the net nucleonic density and shape of the neutron-proton density difference 
are invariant, in the continuous limit, across an isobaric chain. The neutron-proton density 
difference, in particular, is expressible in terms of the inverse operator S for the symmetric 
matter. In this context, we have introduced an isovector density, as a scaled neutron-proton 
density difference, and a counterpart to the net density, also termed isoscalar density. The 
neutron and proton densities can be expressed as a combination of those nearly invariant 
isoscalar and isovector densities. When the local approximation for S holds, the isovector 
density is proportional to the isoscalar density divided by the value of symmetry energy for 
uniform matter. A generalized symmetry coefficient can be introduced for a nuclear system 
and can be expressed in terms of the inverse operator S. The coefficient turns out, further, 
to be inversely proportional to the volume integral of the isovector density. Contributions to 
the net capacitance of the system, for asymmetry, can be expressed in terms of the volume 
integral of the difference of isovector and isoscalar densities. 
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We have carried out Skyrme-Hartree-Fock calculations of symmetric and asymmetric 
semi-infinite nuclear matter, with a significantly higher accuracy than in the past. The ac- 
curacy is important when trying to extract subtle symmetry effects associated with the nu- 
clear surface. We have calculated symmetric surface and surface-symmetry characteristics 
for nearly all Skyrme parameterizations that have been proposed in the literature. In the 
calculations, we have verified the near-invariance of the isoscalar and isovector densities, 
first inferred within the considerations relying on the Hohenberg-Kohn functional. We have 
found that, up to the Priedel oscillations, the isovector density follows the expectation from 
the local approximation, down to about a quarter of normal density. Within the WKBJ 
approximation, we have shown that the isovector density is expected to follow the local 
approximation within the classically allowed region for single-particle wavef unctions. In the 
far-out forbidden region, on the other hand, the density fall-off is governed by separation 
energy. We have extracted displacements of isovector- relative to isoscalar-density as well 
diffuseness values for the two densities. By integrating differences of the two densities for 
symmetric matter, we have determined values of the surface symmetry coefficient for dif- 
ferent Skyrme parameterizations. Consistently with qualitative expectations, we found the 
displacements and volumc-to-surfacc symmetry-coefficient ratios to be strongly correlated 
with the density dependence of symmetry energy in uniform matter. The faster the sym- 
metry energy drops at subnormal densities, the larger the relative displacement of the two 
densities and the larger the coefficient ratio. 

One exciting possibility, emerging as a consequence of our investigations, is that, due to 
the invariance of two fundamental densities, isoscalar and isovector, features of the sym- 
metry energy could be investigated by studying systematics of proton distributions alone, 
without reference to neutron distributions. Such an investigation requires, however, a care- 
ful separation and/or circumvention of shell, pairing. Coulomb and deformation effects on 
the proton distribution for finite nuclei. The obvious difficulty of such separation may be 
compensated by the rather extensive knowledge of proton densities. 
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